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摘要 : Stackel 引力 势 是 一 类 最 普遍 形式 的 可 分 离 势 。 有 具有 Stáckel 形式 的 星系 是 完全 可 积 系 统 ， 
即 其 中 的 恒星 轨道 都 是 规则 的 ， 其 运动 积分 可 以 解析 求 出 。 运 动 积 分 一 一 特别 是 作用 量 (J， 一 
种 特殊 的 运动 积分 )， 能 简化 恒星 运动 的 描述 ; 研究 星系 中 恒星 的 运动 学 和 动力 学 的 重要 途径 是 
作用 量 空间 ， 比 如 使 用 以 作用 量 作为 参量 的 分 布 函数 f(Jr, La, J:)。 通 过 将 一 般 星 系 的 引力 势 
局 部 地 近似 成 Stáckel 势 ， 人 们 可 以 估算 一 般 星 系 中 恒星 运动 的 作用 量 。 介 绍 天 文学 家 在 拓展 
Stackel 引力 势 的 研究 上 取得 的 重大 进展 : 在 基于 Stackel 势 理论 估算 一 般 星系 中 的 作用 量 (或 运 
动 积分 ) 方面 ， 开 发 出 了 若干 快速 数值 算法 (例如 Stáckel 作弊 法 [Stáckel Fudge]. Stáückel 拟 合 
法 )， 基 于 环 面 映射 的 收敛 性 高 精度 数值 算法 (例如 轨道 积分 拟 合 生成 函数 法 、 和 迭代 环 面 构造 法 )， 
以 及 提出 作为 近似 的 运动 积 的 解析 表达 式 ; 在 分 布 函数 的 构造 方面 ， 基 于 上 述 估算 的 作用 量 或 直 
接 利用 Stáckel 势 的 运动 积分 公式 (主要 是 Ia 表达 式 )， 提 出 了 一 些 分 布 函数 模型 ， JR, 工 > Jz) 
或 f(E, Lz, 有)。 这 些 进 展 使 得 人 们 可 以 对 星系 开展 基于 分 布 函数 的 建 模 。 此 外 ， 还 介绍 近 几 年 
把 这 些 方法 应 用 到 银河 系 观 测 数据 的 若干 代表 性 工作 。 

X 键 词 : Stückel 引力 势 ， 作 用 量 - 角 变 量 估算 ;银河 系 分 布 函数 建 模 
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1 引言 和 理论 基础 


银河 系 研究 的 主要 目标 之 一 ， 是 测量 并 计算 出 银河 系 的 结构 ， 即 各 种 组 成 成 分 (如 恒 
星 核 球 、 恒 星 盘 、 暗 物质 晕 和 气体 等 ) 的 空间 分 布 ， 以 及 它们 的 速度 空间 信息 。 为 了 实现 
这 个 目标 ， 在 观测 上 ， 需 要 高 精度 的 天 体 测量 学 及 测 光 学 数据 (如 Gaia) 和 光谱 数据 (如 
LAMOST); 在 理论 和 计算 上 ， 星 系 动 力学 建 模 是 核心 。 星 系 动 力学 建 模 的 理论 基础 是 恒 
星系 统 的 动力 学 ， 即 研究 大 量 的 点 质量 粒子 在 它们 共同 的 引力 场 中 的 运动 "。 关 于 银河 系 
建 模 ， 尽 管 星系 作为 自 引 力 系 统 ， 没 有 星系 能 处 于 严格 的 热力 学 平衡 态 或 完美 的 动力 学 平 
衡 态 ， 然 而 平衡 态 动 力学 模型 (或 称 稳 态 [steady statej， 或 准 定 态 [quasi-stationary state]), 
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对 于 解释 银河 系 和 河 外 星系 的 观测 至 关 重 要 。 从 实际 情况 看 ， 暗 物质 晕 作 为 星系 引力 势 的 主 
要 提供 者 ， 仅 在 短暂 且 罕 见 的 主 并 合 事件 中 ， 才 远离 动态 平衡 ， 一 般 情况 下 该 动态 平衡 仪 受 
吸 积 或 次 并 合 等 事件 的 小 扰动 。 从 模型 需求 的 角度 看 ， 由 动力 学 测量 数据 推断 出 物质 分 布 ， 
该 模型 要 求 假 设 星系 处 于 (HE) 平衡 态 ， 另 外 ， 平 衡 态 模型 是 最 简单 的 模型 ， 更 复杂 的 情况 
可 以 被 建 模 为 平衡 态 模型 的 “扰动 “变形 ”或 迭 加 ”。 

运动 积分 是 粒子 在 轨道 运动 中 的 不 变量 ， 作 用 量 ， 则 是 可 以 与 角 变 量 形成 正则 坐标 的 一 
类 特殊 的 运动 积分 。 由 于 在 “作用 量 - 角 变 量 ” 坐 标 下 运动 方程 形式 的 简单 性 ， 作 用 量 成 为 
很 有 吸引 力 的 变量 。 并 且 作 用 量 无 论 在 微 扰 论 方法 上 还 是 分 布 函数 的 建 模 上 ， 都 是 十 分 必 
要 的 。 但 是 ， 粒 子 在 一 般 引 力 势 (例如 真实 的 星系 ) 中 的 运动 ， 基 本 上 不 能 得 到 解析 的 运动 
积分 或 作用 量 ; 仪 当 引 力 势 为 可 分 离 势 时 ， 能 解析 地 得 到 。Stickel 势 是 椭 球 坐标 系 下 可 分 
离 的 引力 势 形式 ， 该 引力 势 下 的 运动 积分 可 被 解析 地 计算 ”。 在 应 用 上 ， 通 过 将 一 般 引 力 势 
近似 或 拟 合 为 Stackel 势 ， 人 们 可 估算 一 般 引 力 势 下 的 运动 积分 。 近 十 几 年 来 ， 一 些 研究 者 
基于 Stackel 引 力 势 理论 中 的 解析 式 ， 用 (参见 3.2 节 ) 解析 或 数值 的 (例如 Binney 的 工作 ) 方 
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= 法 来 作 拓展 ， 即 基于 Stickel 势 来 估算 一 般 星系 中 的 作用 量 (或 运动 积分 )。 特 别 是 最 近 几 年 
> HK, Binney 和 他 的 合作 者 们 把 他 们 的 方法 都 实现 为 计算 机 程序 ， 并 以 公开 源码 的 方式 释放 


T ( 详 见 第 2 35), 

Stückel 势 作为 一 个 简单 的 引力 势 形式 ， 早 在 19 thee" 就 已 被 给 出 。 在 星系 动力 学 的 
理论 推导 方面 ， 最 普遍 形式 的 三 轴 Stückel 势 在 30 多 年 前 就 已 经 得 到 了 详尽 的 总 结 ， 包 括 
三 个 孤立 运动 积分 (或 作用 量 ) 的 解析 表达 式 ， 恒 星 轨 道 的 分 类 等 "。 但 直到 近 几 年 ， 有 关 
Stickel 势 的 知识 才 在 实际 的 银河 系 建 模 中 应 用 起 来 。 这 主要 归功 于 上 段 所 述 的 Stáckel 势 拓 
展 工作 ， 即 估算 一 般 引 力 势 的 作用 量 (或 运动 积分 ) Te. 

接 下 来 ， 第 1 章 的 其 他 小 节 分 别 介绍 : 恒星 系统 的 动力 学 及 统计 力学 的 概要 、 运 动 积分 
以 及 Stickel 势 的 相关 基本 背景 知识 " 。 第 2 章 介绍 估算 一 般 引 力 势 的 运动 积分 和 作用 量 的 
方法 : 2.1 节 介 绍 Stickel Fudge fg ^ 7* 7^7. 395 节 介绍 Stäcke 引力 势 拟 合法 ™”， 以 
及 更 早期 的 工作 ，2.3 节 介绍 一 类 环 面 坐标 变换 的 作用 量 估算 方法 ， 基 于 Stückel 玩具 势 模 
型 得 到 一 般 引 力 势 下 的 作 有 用量;，2.4 节 是 我 们 基于 文献 [pa] 所 介绍 的 TACT 程序 对 
各 种 方法 做 的 比较 ;2.5 节 则 是 介绍 一 些 其 他 求 运动 积分 的 方法 和 程序 。 第 3 章 介绍 近 几 年 
来 这 些 方法 在 银河 系数 据 上 的 应 用 (如 文献 pup). 38 4 章 进行 讨论 和 总 结 ， 特 别 是 简介 “ 作 
用 量 - 角 变 量 ” 方 法 在 不 可 积 的 星系 现象 ( 见 4.2 节 ) 和 非 平衡 的 星系 过 程 ( 见 4.3 节 ) 中 应 用 
的 最 新 进展 。 

1. ”恒星 系统 的 理论 基础 

本 文 讨论 的 单 粒子 在 整体 势 的 运动 ， 是 从 引力 多 体 运 动 简化 而 来 的 ， 即 单 粒 子 近似 或 称 

平均 场 近似 "。 三 维 背 景 空间 中 一 个 包含 有 N 个 粒子 的 经 典 系统 ， 可 用 哈密 顿 方程 描述 : 


. 0H ~ BH 
Qi 一 Op; , Pi a Óqi p] 


i21,2,3,.,8N . (1) 


对 于 该 系统 所 有 粒子 的 坐标 和 动量 所 构成 的 全 6N 维 相 空间 (一 般 称 为 工 相 空间 )， 通 过 
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相 空 间 代表 


其 中 ， P 


ie, H}pp = 


天 文学 进展 


点 个 数 的 局 部 守恒 可 推导 得 ， 其 概率 密度 pla, p) 满足 刘 
ae 


Ot t io ,Hlsg- 0 ; 
pldi pit), i = 1,2,3, 
Fa OpOH ƏH Op 

| Ogi Opi oqi Op Di 


dp/dt， 即 其 相 体积 不 变 ; 


这 就 是 刘 维 尔 保 体 积 定理 。 


下 面 用 z = 2;,p = p; (j = 1,2,3) 代表 单 


,3N 是 系 综 意 义 的 概率 分 布 ; 
) xo. 左边 其 实 就 是 拉 格 朗 日 视角 的 随 流 导数 
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维尔 方程 : 


T 


p, H 的 泊 松 插 


E 
y 


个 粒子 的 坐标 和 速度 ， 代表 n 粒子 联合 分 布 函 数 。 对 上 述 刘 


f2(qı, Pı, q2, P2, t) 


个 粒子 的 坐标 和 动量 ; dk, Dk AK 


表 第 
维尔 方程 作 单 粒 子 近似 


= fi(qu Pı, t)fi(qs pa;t 加 ， 可 得 到 玻 尔 效 曼 方 程 四 


o 
EE Vif +F: Vaf = | áp f a0 oOo — alfi- ， (3) 
RP, XE f= f(a pt) 是 单 粒 子 分 布 函数 ， 即 处 于 自由 度 3 的 坐标 和 速度 所 构成 的 6 维 


相 空间 (一 般 称 为 相 空 间 ); v 代表 速度 ， 五 代表 受 力 加 ，c(I) 
二 体 碰 撞 ， 对 于 其 他 作用 


考虑 了 
受 方程 体系 


类 型 (如 


中 所 需要 的 假设 非常 繁杂 ( 


刚性 ) 的 散射 ， 将 可 能 


本 文 对 玻 尔 
第 3, 4 章 )。 
TIEA u AR 


效 曼 方程 仅 作 介绍 而 不 作 分 析 (更 多 参见 文献 [B] 中 的 多 


Vs AR quio ir 
eee UR i niente 


Zi. WIK 


um 
o 


vici bel 


此 外 可 以 验证 ， 刘 维尔 方程 在 DAE oD A ERICH A S NRA, RAR 2K 


ZE [RI RI 3-4 PRU CBE A Ze AS RE 0 


对 玻 尔 兹 曼 方 程 进行 局 域 平衡 假设 ， 碰 撞 


项 将 消失 。 


其 实 ， 在 星系 这 种 具体 的 系统 中 ， 


恒星 的 半 


径 与 星系 尺度 相 比 很 小 ， 其 平均 自由 程 远 大 于 系统 尺度 ， 碰 撞 时 间 远 大 于 穿越 整 


个 星系 的 时 


间 ， 比 宇宙 年 龄 还 大 很 多 量 级 ， 因 此 磁 撞 项 本 应 忽略 不 计 。 把 玻 尔 效 曼 方程 两 边 


乘 上 一 个 碰撞 不 变量 E(x v, t) (如 质量 、 动 量 、 


Sb ELA 
He 里 等 )， 


然后 作 局 域 平均 (由 于 碰撞 前 后 的 


可 分 别 得 到 流体 力学 性 质 


玻 尔 效 曼 方 程 两 边 乘 以 


某 些 对 称 性 , 局 域 平 均 后 碰撞 项 不 做 贡献 ， 见 文献 [G0] 第 7A 章 )， 
的 方程 GE: 流体 力学 过 程 是 局 域 平衡 过 程 )。 例 如 对 于 动量 守恒 ， 
£(x,v,1) = mut = x,y, z2)» =o 量 空 间 积分 并 利用 动量 空间 


o 1 

8i a; iti) + Va (pviv) 一 mimm. 
此 式 实质 上 等 同 于 流体 力学 动量 守恒 方程 : 

(pu) _ Pp 


的 散 度 定 理 (无 穷 远 处 分 


(4) 


(5) 


?这 是 一 种 平均 场 近似 ， 即 将 粒子 间 相 互 作用 当成 总 背景 场 ， 而 不 考虑 粒子 间 相 互 作用 入 


联 ; 体现 在 分 布 函数 上 就 是 每 


E: 


粒子 的 单 粒子 分 布 是 无 关联 的 ， 因 此 乘积 为 它们 的 联合 分 布 。 


“由 玻 尔 兹 曼 于 1872 年 提出 。 从 刘 维 尔 到 玻 尔 兹 曼 方 程 的 约 化 ， 中 间 还 有 一 个 BBGKY 级 联 方程 。 

该 受 力 是 对 于 该 “ 单 粒子 ”而 言 的 ， 施 力 者 为 “其 他 粒子 ”( 实 际 上 是 分 布 函 数 ) 的 作用 (因而 在 这 里 直接 体现 平均 场 近似 
的 思想 )， 以 及 可 能 存在 的 外 场 。 玻 尔 兹 曼 方程 可 以 用 来 描述 非 自 治 的 体系 (自治 体系 是 指 粒 子 所 受 力 全 部 来 自 于 研究 系统 本 身 ， 
如 自 引力 系统 )， 即 通过 建立 多 种 组 分 的 玻 尔 兹 曼 方程 来 描述 。 
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其 中 ， 己 为 压力 张 量 。 流 体力 学 的 质量 、 能 量 守恒 也 可 以 类 似 地 得 到 ( 见 文献 [Ed 第 了 章 )。 
这 样 ， 我 们 就 有 了 “ 体系 统 (力学 描述 )- 刘 维尔 方程 (6N 维 )- 玻 尔 效 曼 方 程 (6 维 )- 流 
体力 学 方程 ( 流 场 )” 的 联系 。 
另 一 方面 ， 从 无 碰撞 玻 尔 兹 曼 方 程 出 发 ， 两 边 乘 以 速度 分 量 ， 再 对 速度 空间 积分 并 利用 
散 度 定理 (与 得 到 式 (M) 的 方法 类 似 )， 我 们 得 到 Jeans 方程 : 


pa Oj DG O(vo;;) 
20» A 6 2 7 E (6) 


RP, v 代表 位 形 空间 数 密 SFE, oi; 为 速度 弥散 张 量 。Jeans 方程 不 仅 通 过 平均 简化 了 方程 ， 
还 将 分 布 函数 与 这 些 可 观测 量 联系 起 来 ， 从 而 对 联系 观测 天 文 的 工作 很 有 帮助 。 Jeans 方程 
还 给 人 们 提供 了 一 种 检测 分 布 函数 正确 性 的 方法 。 当 然 , Jeans 方程 仍然 不 是 闭合 的 。 

在 人/ 相 空 间 中 的 粒子 分 布 函数 f = f(a, v, t) 四， 跟 粒子 视角 的 单 粒子 近似 下 的 运动 

(x(t), v(t)) 有 着 密切 的 联系 。 这 种 单 粒子 分 布 函数 ， 实 际 上 假设 体系 是 由 无 穷 多 的 点 粒子 

构成 ， 某 种 程度 上 与 场 论 的 意义 比较 像 ( 场 论 就 是 处 理 无 穷 多 自由 度 )。 从 场 的 角度 来 理解 ， 

这 是 把 所 有 粒子 之 间 的 相互 作用 当做 背景 场 ， 然 后 对 于 任 一 粒子 ， 只 须 研究 它 在 这 个 场 中 

的 运动 ， 不 必 考 虑 其 他 粒子 对 它 的 作用 ; 此 即 前 述 平均 场 的 思想 。 但 是 ， 实 际 上 想 详细 了 解 

星系 结构 与 演化 (如 分 布 函数 )， 直 接 从 自 引力 系统 的 统计 物理 方程 出 发 去 求解 分 布 函数 并 

不 容易 ， 而 且 运用 最 大 入 原理 之 类 的 方法 也 是 不 够 的 (由 于 星系 作为 长 程 相互 作用 系统 的 特 

TE, WB 4 章 )。 为 得 到 分 布 函数 ， 一 种 方法 是 N 体 模拟 ， 基 于 模拟 的 结果 作 统 计 。 另 一 种 

方法 是 根据 粒子 在 引力 势 的 运动 特性 来 构建 分 布 函数 模型 ”， 以 及 其 他 理论 手段 ( 详 见 文献 

R): 在 这 种 方法 中 ， 对 运动 积分 的 研究 是 个 基础 的 环节 。 下 面 ， 我 们 主要 从 哈密 顿 力学 以 

及 单 粒子 近似 下 的 恒星 轨道 运动 的 角度 出 发 ， 介 绍 恒 星 在 静态 引力 势 下 的 运动 积分 的 相关 

内 容 。 

1.2 ”运动 积分 和 作用 量 - 角 变量 体系 简介 

1.2.1 哈密 顿 -雅克 比方 程 

力学 系统 在 相 空间 的 正则 坐标 (q, p) 的 演化 ， 可 以 看 做 是 无 穷 多 个 连续 以 哈密 顿 量 为 

生成 函数 的 无 穷 小 正则 变换 的 结果 。 在 正则 变换 中 令 新 的 哈密 顿 量 为 0， 并 选取 第 二 类 生 


O 成 函数 FA EA ET IE Ais Q, P 均 为 常数 ， 即 — o, P = B， 可 得 正则 变换 的 
方程 : 
Os Os 


被 称 为 哈密 顿 -雅克 比方 程 ( 哈 雅 方程 )， 其 解 的 形式 为 5 = Sla, 6) 就 是 生成 函数 Fl. Hb, 


Ons 


p — mv (m 为 单 粒子 的 质量 )， 而 且 在 星系 动力 学 的 粗 燃 建 模 中 常常 将 测试 天 体 考虑 为 无 穷 多 同 质粒 子 的 分 布 ， 且 
为 方便 一 般 假设 测 试 粒子 的 质量 为 1， 于 是 下 面 我 们 这 里 的 习惯 ， 以 下 讨论 的 了 都 是 关于 æ, v 的 分 布 函数 。 

”第 二 类 生成 函数 是 指 : 在 正则 变换 中 ， 取 独立 的 广义 坐标 变量 为 老 的 坐标 和 新 的 动量 ， 亦 即 生成 函数 为 老 坐标 、 新 动量 以 
及 时 间 的 函数 。 
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第 二 类 正则 变换 中 哈密 顿 函 数 之 间 的 微分 关系 表示 为 : 
Os Os 


= p » = 1,2,3,---, 3 8 
da’ OB? ^ É (8) 


Di 


1.2.2 运动 积分 
如 果 一 个 力学 量 不 显 含 时 间 而 且 保 持 不 变 ，h(gq,p) = const， 就 称 之 为 系统 的 运动 积 
分 ， 即 运动 积分 是 粒子 在 轨道 运动 中 的 不 变量 。 寻 找 系统 的 这 种 不 变量 ， 具 有 重要 意义 。 对 


于 一 个 自治 哈密 顿 系统 ， 哈 密 顿 量 不 显 含 时 间 。 此 时 哈密 顿 量 本 身 就 是 一 个 运动 积分 ， 其 积 
分 常数 马 等 于 系统 总 能 量 。 利 用 公式 0) 和 图 ， 哈 雅 方程 的 解 可 写成 时 间 项 和 广义 坐标 项 


分 离 的 形式 ， 表 示 为 : 

S(q. B,t) = W(q,B) - ht , (9) 
式 中 的 W 称 为 哈密 顿 特征 函数 。 对 于 自由 度 ”的 自 洽 哈 密 顿 系统 ， 最 多 有 7 个 两 两 独立 且 
对 合 四 的 运动 积分 ， 因 此 在 上 式 的 积分 常数 = {PB1,…, Bw} Ahk (n+1) 个 常数 中 ， 至 少 
AEI, BME h = h(61,..., Bn) 相 联系 ， 我 们 取 6, = ho 对 式 (HI) 两 边 微分 得 到 两 
个 微分 关系 ， 以 及 生成 函数 的 微分 关系 ， 即 式 (B): 再 将 W 代入 5 中 ， 可 得 到 关系 W 与 新 
旧 相 空间 坐标 的 微分 关系 : 


| OW a, = OW 
— ôq "o 08i ^ 


现在 我 们 以 哈密 顿 特征 函数 W 作为 生成 函数 作 正则 变换 (gq, p) 一 (Q, P), EF, 
| _OW(q,B)  9W(q B) 


Di (212938 (10) 


OP; Oi (11) 
P; —,, 4=1,2,3,---,n 


此 时 的 新 哈密 顿 量 仍 是 K = A, Jf ABSA BR SP) Mori Aa) EER, 1% 
正则 变换 是 作为 后 面 讨论 运动 积分 坐标 的 基础 。 
更 进一步 ， 对 于 自治 (完全 ) 可 分 离 的 生成 函数 WW， 可 写 为 : 


W = 2 Wi(qi, B) ， (12) 


则 相应 地 可 以 分 离 出 个 哈 雅 方程 (不 对 i 求 和 )， 可 表示 为 : 


OW; 
m (a. FB) = 6. (13) 


哈 雅 方程 方法 正 可 方便 可 分 离 系统 的 求解 。 
MEADE, SAAR n 的 自治 哈密 顿 系统 为 (完全 ) 可 积 的 ， 则 它 存 在 nw 个 互相 对 合 的 


运动 积分 (i = 1,2,3, ,n) (我 们 用 ;代表 已 经 从 E, 8 中 选 好 的 运动 积分 ); 此 即 ( 刘 维 


?如 果 两 个 运动 积分 的 泊 松 括号 为 0， 那么 它们 是 对 合 的 。 


C— 
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AK) 可 积 性 定理 的 内 容 。 另 一 方面 ， 根 据 刘 维尔 保 体积 定理 ， 以 这 些 运动 积分 为 某 一 组 常数 
时 的 坐标 集合 ，M = {a]l p) = 8 ， 它 在 以 五 为 哈密 顿 函 数 的 相 流 下 保持 不 变 ， 并 
HARF n 维 环 面 T" ( 知 不 存在 共振 )， 因 此 将 (完全 ) 可 积 系统 的 轨道 转化 为 一 个 n 维 准 周 
期 运动 ”。 
关于 运动 积分 ， 值 得 注意 两 点 : 第 一 ， 在 星系 动力 学 语 境 中 ， 运 动 积 分 和 运动 常量 是 作 
了 严格 区 分 的 两 个 概念 ; 第 一， 关于 运动 积分 还 有 更 多 细 分 的 概念 对 于 恒星 在 星系 引力 
势 中 的 运动 ， 给 定 引力 场 ， 运 动 常 量 的 定义 是 : 关于 相 空 间 坐 标 和 时 间 的 、 沿 恒星 轨道 不 
变 的 任意 函数 ， 即 C( a(t), vlt), t) = const。 运 动 积 分 也 是 在 恒星 轨道 中 保持 常数 的 量 (HU 
运动 常量 )， 但 它 进一步 要 求 是 不 显 含 时 间 的 函数 ，T( a(t), v(t) ) = const。 运 动 积 分 可 以 分 
为 孤立 积分 和 非 孤立 积分 两 类 ”"。 孤 立 积分 能 影响 恒星 轨道 在 相 空 间 的 分 布 ， 即 相 空 间 轨迹 
(phase trajectory) 所 占据 空间 的 维度 ; 对 于 可 积 系统 ， 每 增加 一 个 孤立 积分 ， 就 把 相 空 间 
轨迹 所 占据 流 形 ( 即 相 空间 的 子 空间 ) 的 维度 降低 一 维 。 非 孤立 积分 则 对 于 恒星 轨道 在 相 空 
间 的 分 布 没 有 影响 ， 因 而 没有 实用 价值 。 
= 另外 ， 在 完全 可 积 的 或 部 分 可 积 的 系统 中 ， 请 注意 对 合 积 分 与 孤立 积分 这 两 个 概念 的 
n 区 别 。 对 于 自由 度 为 双 的 自治 哈密 顿 系统 ， 相 空间 为 2n 维 ， 最 多 有 2n — 1 个 独立 的 ( 即 线 
性 无 关 的 ) 运动 积分 ， 其 中 孤立 运动 积分 最 多 可 以 是 2n — 1 个 ;而 对 合 是 与 可 积 性 有 关 的 
概念 ， 两 两 对 合 的 运动 积分 不 超过 ? 个 ( 当 对 合 积 分 的 个 数 为 n 时 ， 该 系统 就 是 刘 维 尔 当 年 
所 定义 的 完全 可 积 的 )。 因 此 ， 若 有 一 组 运动 积分 是 互相 对 合 的， 那么 它们 肯定 是 孤立 积分 ; 
但 反之 不 然 。 例 如 ， 对 于 球 对 称 引 力 势 中 的 恒星 轨道 人 
孤立 积分 可 以 是 4 个 或 者 5 个 ( 仅 当 轨道 为 开 普 勒 型 时 )， 但 相互 对 合 的 积分 的 个 数 是 3 。 
对 于 星系 动力 学 中 常见 的 轴 对 称 引 力 势 系统 (例如 盘 星 系 )， 一 般 有 3 个 孤立 运动 积分 ， 能 量 
m d EUST RI. Ud dm I. BS CAN 
于 三 轴 引 力 势 系统 ， 一 般 有 3 个 孤立 运动 积分 ， 能 量 E WD h, Is, fH Io, Is 往往 是 没有 解析 
多 式 的 。 过 去 ， 前 人 有 很 多 工作 就 是 在 找 第 2、 第 3 运动 积分 ， 以 及 在 研究 稳 态 系统 的 恒星 
分 布 时 所 用 的 、 以 此 为 自 变 量 的 分 布 函数 f(T) (这 也 是 本 综述 文章 的 中 心 内 容 )。 
1.2.3 作用 量 - 角 变量 体系 
尽管 运动 积分 工 与 其 广义 坐标 ( 记 为 6， 角 变 量 ) 已 经 将 运动 化 为 准 周期 运动 ， 但 是 运 
动 积分 一 般 不 是 正则 坐标 ， 即 它 与 其 广义 坐标 不 满足 哈密 顿 共 斩 ， 运 动 积 分 经 过 变换 后 能 
够 满足 共 罗 的 叫做 作用 量 JJ。 下面 我 们 对 角 变 量 和 作用 量 的 要 求 作 一 下 简介 (可 见 文献 加 
第 3.5.1 节 )。 
如 前 所 述 ， 我 们 假定 存在 从 了 工 到 J 的 可 逆 的 变换 (det(0J/0I) z 0), WEHE J 满足 : 
(1) 运动 积分 ，(2) 正则 坐标 。 于 是 哈密 顿 方程 四 ) 可 用 ， 关 于 J 的 部 分 有 : 


， OH 
TN aay, 14 
Ji 0, 0 (14) 
于 是 哈密 顿 量 只 是 作用 量 的 函数 ， 表 示 为 : 
. A 


~ OU; 
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RP, QJ) 是 角 频 率 。 对 于 束缚 轨道 ， 角 变量 9 和 运动 积分 所 描述 的 就 是 一 个 n 维 准 周期 
运动 ， 即 前 述 的 同 胚 于 AEA HT? 的 流 形 Mi M 上 某 一 维 的 任 一 闭合 环 路 jy; (GEFA 
Fel 51)， 由 茶 一 组 (04, Ji) 所 代表 ， 并 且 要 求人 沿 的 环 路 积分 为 2rr。 

于 是 我 们 可 以 将 相 空 间 轨道 展开 成 传 里 叶 级 数 : 


=) Asie uy X eur, (16) 


角 频 率 之 间 的 比例 关系 丰富 了 体系 周期 类 型 和 遍历 性 的 复杂 性 ( 见 第 4 章 )。 天 体 动力 系统 的 
周期 特性 一 般 来 说 非常 强 ， 这 也 是 作用 量 - 角 变 量 体 系 在 天 体力 学 中 经 久 不 衰 的 重要 原因 。 
定义 作用 量 JI 为 y, 所 围 区 域 的 相 体 积 ”"， 则 有 : 


1 1 1 
2n / / "d dp 2n / LM of yi j dd 2n / I... of ^i i = 
1 
一 j d6;dJ; -4 d0idJi | = Ji- J; , 

T V dud 


其 中 第 2 个 等 式 是 由 于 正则 变换 的 相 体积 不 变性 ， 第 4 个 等 式 是 转化 为 边界 积分 ， 并 且 刨 除 
积分 奇 点 的 邻 域 。 将 奇 点 置 于 作用 量 - 角 变 量 坐 标的 中 心 ，Jf = 0， 此 时 J] 等 于 作用 量 的 不 
变量 。 于 是 ， 我 们 由 正则 变换 的 庞 加 莱 不 变量 的 不 变性 ， 可 得 : 


—— id0 = 一 , 1 
J TEE x. f, nda (18) 


其 中 ,人 是 正则 变换 的 矩阵 ， 指 在 y 对 应 的 相 空间 所 在 区 域 。 也 由 此 可 见 ， 某 个 作用 量 的 
值 .代表 子 环 y: 所 围 的 面积 ( 当 对 qu, p 选取 合适 的 标 度 )， 而 0; 代表 转 过 的 角度 。 

KEE ER 以 及 式 (D 的 第 1 式 ， 易 得 ps = SPREE at ec 
系统 ，p 只 依赖 于 相应 的 gi 而 独立 于 其 他 广义 坐标 。 这 样 一 来 ， 作 用 量 的 计算 可 简化 为 
J= 5 $ pidge 并 且 Wi Br std 8 运动 一 局 而 其 他 不 变 ， 回 到 初始 状态 ， 此 过 程 中 哈 
密 顿 特征 函数 W, 的 变化 量 AW; = pidgq;， 正 好 是 作用 量 的 2 倍 ， 这 就 是 作用 量 的 意义 。 
此 外 ， 我 们 可 用 式 (5) 计算 角 频 率 ， 用 式 (ID) 中 第 1 行 的 哈密 顿 特征 函数 计算 角 变量 
的 值 。 数 值 计算 上 见 附录 A. 

于 是 ， 从 一 般 的 相 空间 坐标 到 作用 量 的 关系 为 ， 相 空间 坐标 (q, p) -> 运动 积分 工 和 角 
BERO 一 正则 环 面 坐标 (J, 0). 
12.4 ”基于 守恒 量 的 分 布防 数 

关于 单 粒 子 分 布 函数 f (oo, v) (此 时 为 3 个 自由 度 的 ， 由 坐标 、 速 度 所 构成 的 6 维 p 相 空 
间 )， 对 于 QE) 稳 态 体系 ( 指 9f/8t = 0)， 由 无 碰 挤 玻 尔 兹 曼 方程 可 推导 ， 分 布 函 数 可 写成 
运动 积分 的 函数 f — f(D). HERD Jeans 定理 ， 而 对 于 规则 轨道 具有 非 谐振 频率 的 准 稳 态 恒星 
系统 的 分 布 函数 ， 可 写成 3 个 作用 量 的 函数 f = fJ) HEIR Jeans 定理 。 因 此 人 们 能 
加 不 元 余 、 深 刻 地 理解 恒星 的 分 布 。 


(17) 


imd 
T 
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早先 前 人 所 用 的 分 布 函数 模型 ， 常 常 是 FUE.) 的 函数 形式 ， 即 把 能 量 作为 分 布 函数 
的 参数 (参见 E23 节 )。 不 幸 的 是 ， 对 于 实际 的 建 模 工作 而 言 ， 这 是 一 个 陷阱 ， 从 而 阻 滞 了 基 
于 运动 积分 或 作用 量 建 模 这 项 事业 的 进展 ”。 对 于 实际 星系 的 动力 学 建 模 ， 不 能 将 E UN 
分 布 函数 的 参数 列表 中 。 原 因 很 简单 : 实际 星系 是 一 个 多 组 分 系统 ， 动 力学 模型 中 至 少 要 包 
括 暗 物质 晕 和 恒星 组 分 。 轨 道 的 能 量 不 是 局 域 量 ， 如 果 五 出 现在 参数 列表 中 ， 那 么 该 轨道 
周围 的 相 空 间 密 度 将 以 不 受 控 制 的 方式 变化 。 例 如 ， 最 简单 的 轨道 是 恒星 位 于 银河 系 中 心 的 
规则 轨道 ， 当 添加 暗物质 党 组 分 后 ， 中 心 引力 势 会 变 深 ,位 于 中 心 的 恒星 的 能 量 会 减少 ， 而 
该 轨道 保持 不 变 。 但 是 ， 只 有 当 模 型 完全 组 装 好 并 且 确 定 了 其 势能 之 后 ， 我 们 才能 知道 该 轨 
道 的 能 量 Eo = 8(0)。 因 此 我 们 无 法 在 建 模 工 作 一 开始 ， 就 指定 恒星 的 f(E,---). Taba 
函数 只 是 作用 量 的 函数 时 ， 每 个 轨道 周围 的 相 空间 密度 则 不 会 改变 ”。 

总 之 ， 相 比 于 一 般 运 动 积分 ， 作 用 量 - 角 变 量 体 系 的 优势 是 : 作用 量 是 一 种 特殊 的 运动 
职 分 ， 而 且 作 为 正则 坐标 ， 作 用 量 有 共 斩 坐 标 即 角 变量 ， 作 用 量 体 系 适用 于 使 用 微 扰 方 法 ， 
也 正 适 合 描述 周期 性 强 的 天 体 运 动 ， 作 用 量 是 绝热 不 变量 ， 可 用 于 缓慢 变化 的 引力 势 ， 以 及 


= 邻近 势能 中 的 轨道 识别 ， 此 外 如 上 段 所 述 ， 使 用 以 作用 量 为 参数 的 分 布 函数 来 为 星系 (多 组 
> 分 系统 ) 建 模 ， 才 使 得 基于 分 布 函数 的 建 模 成 为 可 能 。 


1.8 ”Stickel 势 简介 

Stückel 势 在 椭 球 坐标 下 可 分 离 @。 三 轴 Stackel 势 是 一 类 最 普遍 形式 的 可 分 离 势 。 粒 子 

在 一 般 Stáckel 势 中 的 运动 ， 能 通过 代数 手段 得 到 解析 的 运动 积分 ， 即 是 完全 可 积 的 。 本 节 
介绍 三 轴 Stickel 势 的 基本 内 容 ”， 在 这 之 前 先 回顾 下 椭 球 坐标 系 。 

作 坐 标 变换 (x,y,z) 一 (和 ,14,v)， 并 令 

ja y? 52 

Tta T+B 7 十 了 yy H 

v 即 构成 共 焦 椭 球 坐标 。 令 常数 a = a, B= —b?, y = -c?, FF 

—y € v € —B € p € —a € M 47 变化 时 ， 半 轴 长 变化 但 焦距 不 变 ， 形 成 一 系列 


(19) 


的 共 焦 二 次 曲面 。 该 三 维 坐 标 系 的 等 坐标 面 如 图 1 所 示 ， 如 式 (TH) Fer =r. AF 
(o à+a >20, 和 十 6 之 0, 和 +Y 关 0， 此 时 每 个 等 入 曲面 代表 一 个 椭 球 面 (图 1 中 的 红色 面 ); 类 
= fU, HF uta<0,p+B>0, n+ 7>0, Sp HRA OH (黄色 ); 等 > 曲面 
© 代表 双 叶 双 曲 面 (绿色 )。 和 卡尔 坐标 系 、 球 坐标 系 等 ， 可 看 作 共 焦 椭 球 坐标 系 的 极限 情形 。 
u 共 焦 椭 球 坐 标 系 也 是 正 交 坐标 系 ， 其 线 元 ds? = Yo, Pdr, HEMRA: 
P= (A= p) =v) 
40+ o)0  B)0 3) 
(p — v)(u — A) 
P= 20 
Au + a)(ut B)(u+7) en 
P= (v - Av - p) 


” 4(v-a)(v + Bv +7) 


?Stückel "4tEff] Ll EUEWI, H = dp? + (a) 的 哈密 顿 -雅克 比方 程 分 离 的 唯一 坐标 系 是 共 焦 椭 球 坐标 系 (参考 文献 [E] 第 
3.5 节 )。 
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ik: 椭 球 面 (红色 ) 是 等 和 面 ， 单 叶 双 曲面 (黄色 ) 是 等 y 面 ， 双 叶 双 曲面 (绿色 ) AEE vq. AB ba T EK 
度 比 例 。 


1 ， 共 焦 椭 球 坐 标 系 的 等 坐标 面 示意 图 


此 时 Stackel 势 5s(X, u,v) 可 写成 下 面 形式 : 
CA) nu) El) 
Pop P 

F(A) F(u) F(v) 
(A-u)A-v) (u-v)u-AX) (v-uv-nu) 
FEA, Cn, 6 分 别 是 相应 椭 球 坐标 和 ,1,v 的 任意 函数 ( 见 文献 团 中 的 式 (6)) B. + F(A) 
A(X + a)(A+ B) - 3)C(). F (u) K F (v) KIERRE, BBA Stickel 势 可 以 写成 3 个 分 量 
分 离 的 形式 。 

将 上 面 的 Stáckel 势 代 入 哈密 顿 量 


2 2 
Dx Pu Dy 
一 | | FOs(A,n,v) , 22 
2P2  2P2 2P? ss Hy) "A 


Bs 


, (21) 


E 


?不 同文 献 对 引力 势 所 采用 的 符号 有 时 候 不 相同 (V 0 或 其 他 )， 相 应 地 有 些 文献 中 哈密 顿 量 记 为 H = ip? 十 了， 有 些 记 
AH = 3p? 十 $B; 类 似 地 ， 对 其 分 离 函 数 F (r) 前 面 的 正 负 符 号 的 规定 ， 也 可 能 有 所 出 入 。 本 文采 用 如 式 (QO) 所 示 的 方式 。 

?为 了 书写 方便 ， 本 文 后 面 统一 用 F(r) 代 指 分 离 函数 F(A) 等 3 个 分 离 函数 ， 而 实际 上 这 3 个 函数 可 以 采取 不 同 的 具体 
形式 。 
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其 中 ，p、 = PPA p, Rp, 依 此 类 推 。 再 代入 哈密 顿 - 雅克 比方 程 瑟 = E, 以 及 有 
p. 二 0W/97， 其 中 W 为 哈密 顿 特征 函数 。 解 哈密 顿 -雅克 比方 程 ， 令 两 边 乘 以 (和 一) (4 一 
v)(v 一 入 )， 可 得 


(v = n) * a)(A+ B) t ) 


( — V) B(u + o) (a + B)(u 3) ) F(u) — pE] (23) 


(u — A)[2(v +a) 4- B)(v 4- y) ( 5) F(v) -CWE| «0. 


juni 


HFR (Er) 的 特殊 形式 以 及 W 的 可 分 离 性 ， 可 写 为 W = Z,a Wi(7)。 于 是 式 (3) 中 
括号 内 的 部 分 只 是 各 自 7 坐标 的 函数 ， 设 之 为 U;(7) 并 蔡 换 ， 可 得 : 


(v = WDAN+tTO- Uy) + (u = A)U w) = 0 . (24) 


将 上 式 两 边 对 每 个 坐标 r = {A uv} 作 两 次 微分 ， 即 将 on 作 用 于 式 (Za) 两 边 。 然 
后 我 们 从 其 中 的 对 同 坐 标 二 次 微分 式 中 可 得 oUr) = 0， 即 3 个 U.(r) 是 线性 函数 
U.(r) = jer 一 kr; 利用 不 同化 标的 二 次 微分 式 ， 可 得 3 个 U (T) 的 线性 系数 相等 ，j; = j; 
K U (T) RAR a), WR OA, u,v} 取 值 的 任意 性 (或 者 再 对 每 个 坐标 T = A u, v 作 微 分 )， 
可 得 = 相等 ，k; — k, BB 7 大 就 是 可 分 离 方程 的 分 离 常数 。 因 此 3 个 U(r) 只 能 是 相同 
的 线性 函数 


UAT) S7F = hs (25) 


HF U, (7) 即 式 (E3) 的 中 括号 内 的 部 分 ， 则 可 得 : 


juni 


à (Amy. TE-cTjck-F(r) 
2( 


fr Or THal(r+B)(r +7) ’ (26) 


此 即 Stáckel 势 下 动量 和 运动 积分 J,K 关系 式 。 注 意 ， 常 数 j, 就 是 除 哈 密 顿 量 之 外 两 个 
运动 积分 J, K 的 值 (类 似 于 总 能 量 和 哈密 顿 量 H 的 关系 )， 后 面 本 文 我 们 不 再 刻意 区 分 ; 
类 似 的 还 有 do, da 之 于 D, Da. 

EHE, H, J, K 可 写 为 : 


H=X+Y+Z 
J=(u+v)X -(v-3MY -(A-u)Z , (27) 
K = uv X 十 ZAY + AuZ 


| (ERHTIJ 


ChinaXiv 合 作 期 刊 


154 天 文学 进展 39 3$ 
其 中 ， 2 

x= Dx F(A) 
2P? ^ (X= p)(—2) 
py F(p) 

Y = H 28 
2P?  (u—v)(u- 2) (28) 
r F(v) 


此 外 ， 还 有 更 常用 的 另 一 套 运 动 积分 ， 两 者 的 关系 是 : 


o?H +aJ+K 
h= 
"p (29) 
yH+y+k 
L3 = ———————— 
y—a 


JC BUR AAP BO BSE” Da 1s。 相关 推导 见 文献 团 。 类 似 
Hh, Hip. = 9W/97， 得 到 动量 的 表达 式 为 
^ quc PONE renis Ure COH PUR) 
2(r + o)(r + B)(T +7) l 
上 式 正 是 用 来 计算 Stickel 势 下 的 运动 积分 ， 也 是 估算 一 般 势 的 运动 积分 时 主要 用 到 的 动量 
表达 式 。 由 作用 量 的 定义 ( 即 式 (OR), FIA: 


2 fF 
J, = Ọ pdr = =| \pr|dr . (31) 


积分 即 得 作用 量 ， 通 过 哈密 顿 特 征 函 数 则 可 求 出 角 变 量 9。 附 录 A 介绍 了 如 何 具体 数值 求 
解 作用 量 和 角 频 率 ( 见 文献 中 221). 
此 外 ，Stickel 势 还 有 一 个 优势 ， 当 星系 处 于 准 稳 态 时 ， 恒 星 分 布 函数 的 非 对 角 二 阶 速 
度 矩 消失 ， 从 而 仅 剩 的 关于 三 个 对 角 二 阶 矩 的 3 个 Jeans 方程 变 成 了 闭合 的 ; van de Ven 等 
人 四 用 微分 方程 的 方法 给 出 了 Stückel 势 下 Jeans 方程 的 速度 弥散 通 解 。 


2 基于 Stickel 势 近似 的 作用 量 估算 方法 


前 面 介绍 了 作用 量 - 角 变量 的 意义 以 及 Stackel 势 下 运动 积分 的 解析 求解 ， 然 而 绝 大 多 
数 实际 引力 势 下 并 没有 解析 的 运动 积分 。 接 下 来 本 章 介 绍 几 种 估算 一 般 引 力 势 系统 的 作用 
量 - 角 变量 的 估算 方法 ， 这 些 方法 合作 用量 - 角 变 量 体 系 得 以 应 用 于 一 般 引 力 势 系统 。 由 于 借 
己 知 的 某 种 具体 形式 的 Stáckel 势 作为 玩具 模型 ， 因 而 这 些 方法 其 实 列 含 了 微 扰 法 的 思 
想 。 具 体 来 说 ， 这 些 方法 的 过 程 是 : 给 定 某 种 已 知 的 具体 形式 的 Stickel 势 作 为 玩具 模型 势 ， 
并 在 一 定 程度 上 用 所 画 的 “真实 ”轨道 上 的 采样 点 数据 来 拟 合 相 应 的 待定 模型 引力 势 (或 者 
生成 函数 里 的 参数 )， 然 后 用 这 个 计算 的 势 作为 “真实 Stückel 势 ” 来 估算 运动 积分 ， 从 而 计 
算 作用 量 及 角 频 率 等 所 需要 的 量 。 
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Sanders 和 Binney ”还 提供 了 相应 的 开源 程序 TACTE, 该 程序 用 C++ 所 写 ， 能 
够 用 若干 种 方法 计算 出 作用 量 、 角 频率 和 和 角 变量 ， 并 能 给 出 误差 ， 另外， 还 提供 了 相应 的 
Python 画图 程序 。 关 于 算法 的 更 多 细节 ， 请 参见 文献 U,Zd. WR TACT 程序 的 源码 。 
2.1 Stáckel (FRE 

关于 估算 一 般 引 力 势 系统 的 作用 量 - 角 变量 ， Binney” 提出 了 高 效 的 Stáckel Fudge 77 
法 ， 可 看 作 是 一 种 作 浆 算法 (以 后 称 为 Stáckel (EMA). FARE: 假设 引力 势 可 以 局 部 地 
写成 Stáckel 势 的 形式 ， 那 么 以 一 些 技巧 局 部 地 ( 即 在 恒星 轨道 的 某 一 段 之 内 ) 计算 目标 势 的 
伪 分 离 函 数 (参见 式 (ZT) 的 R(r))， 并 且 假定 其 他 坐标 对 分 离 函 数 误 差 的 交叉 影响 较 小 ， 然 
后 以 此 计算 一 组 新 的 运动 积分 ， 从 而 通过 积分 得 到 作用 量 。 下 面 先 简介 Binney 
首先 提出 的 轴 对 称 情形 的 Stickel 作弊 法 ， 然 后 介绍 三 轴 Stáckel (EME, ， 后 者 是 对 前 者 
的 扩展 。 
2.1.1 轴 对 称 势 的 Stäckel 4E HEE 


LT 


一 般 来 说 ， 轴 对 称 势 下 有 3 个 孤立 运动 积分 ， 即 总 能 量 E. AIEEE LS UE 
v J。 轴 对 称 椭 球 坐标 系 相对 于 三 轴 椭 球 坐 标 系 而 言 ， 一 个 半 轴 长 参数 与 男 外 一 个 相等 ， 如 果 


ra a 二 6， 则 为 长 等 轴 (prolate) 椭 球 坐标 系 (A, v) HPF o 代表 某 个 子午 面 的 方位 角 ， 如 果 
< 十 y= 6， 则 为 扁 等 轴 (oblate) 椭 球 坐标 系 (A, i, x)». rh x 代表 某 个 子午 面 的 方位 角 。 因 此 
在 描述 粒子 在 轴 对 称 势 中 的 运动 时 ， 可 以 消去 一 个 自由 度 。 

基于 等 轴 椭 球 坐 标 系 作 出 新 的 (u,v, p) 坐标 (参考 文献 加 第 3.5.3 7), HAR = 
Asinhusinv, z= 二 和信 coshucosv。 该 坐标 系 下 轴 对 称 Stáckel AWN: 


F,,(u) B F,(v) 


$s(u,v) = : 32 
s( sinh? u + sin? v (32) 


由 类 似 于 式 (BU) 的 关系 ， 需 要 显 式 的 Stáckel 分 离 函 数 也 (7)， 才 能 得 到 用 于 计算 作用 量 的 
动量 pz. 

E HRI AA O(u, v) 写成 Stackel 势 的 形式 ， 则 E (u)— F, (v) 
= 并 且 知 目标 引力 势 接近 Stackel A, M 6E, (u) 对 wv WMA OF, 


x (sinh? u--sin? v)ós(u, v); 
(v 


) 对 v 的 依赖 是 很 小 的 。 


OF,,(u) = (sinh? u + sin? v)d(u, v) — (sinh? uo + sin? v)®(uo, v) zz F,(u) — Fu(uo) , (33) 


同样 ，F,(v) 可 以 重 写 为 及 (rr/2) + 8F, (v). 
这 里 wo 是 一 个 不 影响 作用 量 计算 的 参考 值 ， 其 选取 是 固定 v 而 使 àF,(u) 最 小 值 时 的 
( 见 文献 Ia)。 已 (wo) IF, (7/2) 的 计算 为 : 


Phttps://github.com/jls713/tact 
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2 I? 
I; + F, (uo) © Ia + F,(u) — ôF, (u) = E sinh? u eu = 


Es 
2^?  2^2sginh?u (u) 
2 2 
. p L 
I, + F,(n/2) z Iz + F,(v) - 6F,(v) = —Esin? v + 2> 4 : F, 
3 ( / ) 3 (v) (v) in v 2A2 2A2 sinh? v (v) 


(34) 
2.1.2 三 轴 执 的 Stäckel FH 
本 小 节 主 要 介绍 Sanders 和 Binney 提出 的 算法 ， 其 是 对 Binney ™ 提出 的 Stickel 作 
次 法 的 三 轴 推 广 。 对 于 一 般 轴 对 称 势 系统 ， 除 了 Stickel (EWE, Sanders 还 提供 了 对 一 
般 轴 对 称 势 的 Stáckel 拟 合 的 方法 (下 一 小 节 介 绍 )。 但 是 三 轴 Stáckel 势 比 轴 对 称 的 多 一 维 
参数 ， 对 其 拟 合 所 需 的 具体 拟 合 模型 会 复杂 许多 ， 于 是 Sanders 和 Binney 推广 出 了 三 轴 势 
的 Stackel 作弊 法 。 该 方法 借助 类 似 于 运动 积分 J, K 的 另 一 套 运 动 积 分 (作弊 法 得 到 的 新 的 
不 独立 的 运动 积分 )， 以 及 局 部 近似 的 Stáckel 势 ， 用 动量 和 运动 积分 的 关系 式 得 到 每 一 点 的 
动量 pr， 再 积分 得 到 作用 量 J; 类 似 地 ， 可 求 出 角 变 量 9。 本 方法 以 当前 真实 的 相 空间 点 做 
pus 初始 条 件 (xo, vo)， 来 计算 一 系列 动量 ， 这 样 就 能 不 需要 Stáckel 近似 势 的 分 离 势 Fr) 的 显 
< 式 表达 式 ; 然后 积分 得 到 作用 量 ， 最 小 化 作用 量 的 标准 差 (与 每 个 初始 (xo, vo) 有 关 )， 调 整 
焦距 ， 从 而 再 得 到 合适 的 作用 量 ; 对 于 轴 对 称 势 系统 ， 还 可 做 迭代 以 减 小 误差 。 下 面 介绍 其 
LACER. 
(1) 三 轴 势 Stáckel 作弊 法 与 运动 积分 的 计算 
对 于 一 个 一 般 引 力 势 ， 定 义 三 轴 情 形 的 辅助 函数 : 


XA. p, v) - (A = u)(v > A)9(A, p, v) 
X4 Os H, v) = (u a v)(A = L)P(A, H, v) ’ (35) 
Xv (A, 4, v) (v E A)(u z v)9(A, u,v) 


如 果 引 力 势 就 是 Stáckel 势 ， 那 么 式 B 可 以 写 为 : 
F(u) - FY) , vF(u) - uF(v) 


xalA, Hs v) = F(A) 入 "-— | u-v (36) 
现在 ， 假 设 该 引力 势 是 近似 于 Stückel 形式 的 ， 从 而 由 上 式 可 以 写成 : 
F(T) = x«(A. uv) + TC» + D. . (37) 


当 固定 另外 两 个 酉 球 坐标 时 ， 某 个 r 对 应 的 Cr, D, 可 看 作 常 数 ， 从 而 使 我 们 可 以 分 开 来 独 
立 计算 所 需 的 运动 积分 。 
将 式 E 代入 到 式 (Z6)， 得 到 动量 的 另 一 个 计算 式 


2(r+o)(r+P)(r c pr = TE - TJ, + Kr + Xr(A MY) - (38) 


与 式 (E) 相对 比 ， 就 有 了 这 些 常 数 之 间 的 关系 式 J, = 7 -—C,, K, = k+ D, (j, k 正 是 前 面 
A (2g) 的 运动 积分 )。 这 样 ， 给 定 一 个 初始 相 空 间 点 (zo, vo) 和 一 套 合适 的 椭 球 坐标 系 (AK 
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标 变换 参数 即 (o, 6,Y))， 以 及 引力 势 各 点 的 值 ， 便 可 计算 某 椭 球 坐标 的 相 空 间 点 时 的 运动 


py. 


下 面 介绍 J,, K, 的 计算 。 由 式 (Ga) 得 : 


K, = 2(% + a)(70 + 8)(70 4 qr, TEE + To — X« (Ao, Ho, vo). - (39) 


K, 依赖 于 万 ， 而 其 他 量 可 以 计算 。 当 固定 另外 两 组 椭 球 坐标 和 动量 时 ， 由 于 Stückel 势 的 
分 离 性 ， 可 以 通过 仅 改变 某 个 > 和 p, 的 值 ， 便 可 使 粒子 仍 保持 在 既定 轨道 上 ， 从 而 认为 在 
Stickel 势 下 保持 椭 球 相 空 间 坐 标 不 变 ， 哈 密 顿 量 对 某 个 椭 球 坐标 的 偏 导 为 0 ( 见 文献 四 中 
的 第 3 章 )， 即 ; 


_ OH q0 [m 1 mo ,1 È 3 
0= Cay)" = 33 E] | 2(u—A)P2 2(v—X)P2 DA ay) 
分 别 将 入 和 Ao RAR (BG), JEVEZE, KRA Po UTES 
pu 1p _ [200 + a) (Ao + B) (Ao + $)p3, + XQ: Hos vo)] + (A? — AE — (A= Ao) Ja 
N 2B (A= p)(A-v) (41) 
Bel P(A, Lo, Yo) . 
{xX (ED) 来 计算 式 (uu) 第 一 项 ， 可 见 同 时 可 以 消去 96/6A， 因 此 得 到 Ja 以 及 类 似 的 Ju. Jy 
的 表达 式 : ip. 2 
1 Pig (Ao H 1 py, (Ao Y 
Jy = (Ho vo)E | 2 £ Pi ° 2 pi 2 
_ 1px, (uo — Ao) , 1%, (uto — vo) 
Jy = (Ao + vo) E 4 2 : PE : 2 PE - (42) 
1 p3, (vo — Ao) | 1 pz, (vo — uo) 


T 


Jy — (A E 
(Ao + Lo) t5 P; 2 P2 


于 是 得 到 6 个 估算 的 运动 积分 J., K, 的 表达 式 ， 其 中 只 有 3 个 是 独立 的 ” 。 
(2) 合适 焦距 的 选取 
(3 有 了 这 些 新 的 运动 积分 ， 可 用 式 (B8) 求 得 相应 7 点 的 动量 ， 之 后 积分 得 到 作用 量 。 
= 轴 对 称 情 形 也 类 似 。 然 而 ， 坐 标 变换 的 参数 的 焦距 (A, = VB 一 a, As = V7 二 5， 任 取 
y = -1kpc)， 其 大 小 非常 影响 该 方法 所 估算 的 作用 量 的 方差 (认为 是 误差 )。Sanders 和 
Binney 给 出 了 数值 调整 合适 焦距 的 方法 。 例 如 对 于 短 轴 闭环 轨道 (short axis closed loop 
orbits)， 轨 道 形成 在 z = 0 面 ， 焦 距 为 y = 土 A1 的 椭圆 。 在 这 种 轨道 ， 只 有 作用 量 J, 不 为 
0， 利 用 这 一 点 ， 给 每 个 能 量 值 找 出 一 个 合适 的 焦距 。 
其 步骤 是 ， 先 给 定 一 个 一 般 势 ， 最 大 和 最 小 能 量 值 以 及 期 间 的 能 量 网 格 点 ;在 每 个 网 格 
点 E; dHB—^ ^ MAB ER (0, yi, 0) 以 速度 为 V2(E; — ®(0, yx, 0)) 垂直 出 发 的 半 周 期 的 轨道 ， 记 
录 所 画 轨 道 的 相 空间 点 并 计算 作用 量 ,， 求 出 使 作用 量 的 方差 最 小 的 焦距 A AS 的 计算 
也 可 用 上 述 方法 求 出 。 
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(3) 环 面 迭 代 增 加 精度 
Stackel 作 次 法 的 过 程 复杂 ， 且 误差 较 大 ; 偏差 范围 可 达 百 分 之 几 至 几 十 ， 特 别 是 盒 轨 
道 的 偏差 更 大 ( 见 文献 [zz] 中 的 图 6 一 8)。 但 是 若 把 它 与 环 面 映射 ( 见 区 3 节 ) 方法 相 结 合 ， 
通过 环 面 兴 代 ， 可 使 误差 收敛 。 这 种 迭代 算法 将 在 区 3 可 市 具体 介绍 。 
2.2 Stickel 拟 合法 
2.2.1 一般 轴 对 称 势 的 Stáckel MAA 
这 里 介绍 Sanders" 提出 的 长 等 轴 椭 球 坐 标 系 下 的 一 般 轴 对 称 势 的 Stickel 拟 合法 。 该 
方法 是 对 一 般 引 力 势 进行 Stáckel 形式 的 拟 合 ， 即 采样 轨道 来 拟 合 出 一 个 合适 的 Stáckel 
引力 势 的 分 离 函 数 F(r)， 来 得 到 运动 积分 。Sanders 和 Binney 四 把 Stickel 作 浆 法 和 这 种 
Stackel 拟 合法 归 类 为 “ 非 收 敛 方法 ” 因为 这 两 种 方法 依赖 于 做 一 定 的 近似 ， 而 近似 所 导致 
的 误差 不 能 通过 增加 计算 量 来 不 断 减 小 。 
在 某 个 子午 面 上 ， 运 动 的 位 形 空间 从 柱 坐 标 (R, o, z) 变换 到 扁 球 坐标 (A, p, v) 依赖 于 
R? g? 
wy THQ Li qp 一 
其 中 ， 和 ,了 是 上 式 作 为 关于 r 的 方程 根 ，@? = a, bP = B, c = yY， 类 似 于 三 轴 情 形 下 的 式 
(9)。 此 时 Stäckel 势 的 形式 简化 为 ; 


1, (43) 


S Z0 - Fo) 
$s--— e (44) 


轴 对 称 时 第 三 运动 积分 的 解析 表达 式 为 : 


I2 
2(A + a) 


F(A) XA2(A— wv)? 
入 十 了 8(Ata)(A+7)? 


(45) 


动量 和 运动 积分 的 关系 表示 为 : 


; 1 Ls F(r) Ts 
Pr 2(7 +a) E aa) TFT >] : (46) 


o 下 面具 体 介 绍 对 一 般 轴 对 称 势 的 Stáckel 拟 合 和 环 面 坐标 建立 的 方法 以 及 每 步 所 需要 的 


(1) 给 定 初始 的 (z,v)， 以 及 真实 的 势 (轨道 积分 中 所 感受 到 的 势 )， 并 在 真实 的 势 中 计 
算 实际 的 运动 积分 E, Lz 

(2) 焦距 的 选取 。 作 上 述 变 换 到 椭 球 坐标 。 然 而 ， 变 换 的 参数 ， 即 焦距 A = a? — c? ( 因 
为 a 二 2， 只 要 一 个 2 一 局 便 可 确定 椭 球 坐标 系 ， 可 取 c = 1)， 非 常 影响 作用 量 的 标准 差 ， 
因此 通过 一 定 的 手续 来 估算 一 个 更 好 的 焦距 A. Sanders" 采用 的 估计 方法 是 ， 由 于 Stückel 
势 的 分 离 性 ， 可 得 


82 


3x3; (^ — 9s) - 0 . (47) 
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由 上 式 以 及 两 种 坐标 之 间 的 变换 关系 ， 可 推 得 ” : 


OFs Obs | p, 4 =] /Ss 


2 2 p 2 
dew eder est E 0m dz OROz ` 


(48) 


既然 假设 所 估算 的 局 部 势 D 就 是 该 Stickel WA 5s， 于 是 就 用 该 式 来 估计 a? 的 值 ， 其 中 的 
偏 导 用 已 知 的 真实 更 来 数值 计算 梯度 而 得 ; 每 个 a? 值 ， 需 要 在 轨道 中 每 个 点 的 相应 值 来 计 
算 ， 最 后 做 平均 ， 这 些 点 在 下 一 步 的 轨道 采样 中 获得 。 

(3) 轨道 采样 以 及 运动 区 域 的 估算 。 进 行 轨道 积分 ， 即 在 (z,v) 坐标 下 通过 每 一 点 的 
势 和 时 间 积 分 来 算 下 一 步 的 坐标 (an^, v)， 即 相当 于 数值 法 画 轨道 。Sanders 和 Binney 的 
TACT 程序 中 用 守恒 格式 的 数值 方法 画 出 轨道 ， 并 且 每 个 轨道 周期 离散 为 300 个 点 。 

在 画 轨道 时 : 1) 在 前 几 周期 的 一 些 点 中 ， 隔 5 个 点 用 式 (GS) 估算 a? (c? OR 1) 并 平均 ， 
最 后 得 到 那个 最 合适 的 扁 球 坐标 ; 2) 通过 前 后 Or/Ot 异 号 等 条 件 (其 中 每 次 需要 zz 到 7 的 坐 
标 变换 ) 来 判断 是 否 达到 轨道 的 7 的 端点 ， 找 到 后 就 记录 ， 得 到 和 vuv = c)» Matt 
到 轨道 在 椭 球 坐标 的 所 在 区 域 。 
(4) 样 条 插值 拟 合 Stückel 势 。 关 于 引力 势 拟 合 ， 这 里 假设 局 部 是 Stickel 的 ， 于 是 
先 在 椭 球 坐标 下 若干 点 计算 出 相应 的 “假设 是 Stackel 形式 的 所 求 势 ” 的 值 ， 之 后 对 这 
些 点 的 势 的 值 做 三 次 样 条 插值 得 到 任意 点 的 值 。Stickel 辅助 函数 x 在 轴 称 下 简化 为 
x= (A, v)B(A, v). 

假如 所 拟 合 的 势 5 是 Stackel WIN, BA A, vB, v) = F(A) 一 了 (v)， 正 是 据 此 得 到 
若干 插值 点 局 部 的 假设 是 Stückel 形式 的 所 求 引力 势 的 值 。 关 于 F(r) 的 选取 ， 其 要 求 也 正 
是 使 其 与 Stackel 势 的 偏差 


Bip] = [ax f° avapyN (A0 0 = FQ) + FW)? (49) 


最 小 。 归 一 化 的 权重 函数 形式 为 : 


ACA) 24A (A AAT, N) = (g -rv )!. (50) 
这 个 权重 函数 是 提前 选 的 ， 相 关内 容 可 见 文 献 四 , 团 。 于 是 拟 合 中 所 选用 的 Stáckel 势 的 具体 
PRN: : : 
PIS XO) Xy Po) =e) Ss (51) 
其 中 ， 
Vy 入 十 
A= f XANGO), xe) = f° axo 40) . 
Me ies - (52) 
z= / a f dvx(A, v)A(A) N(v) . 
这 是 因为 在 这 种 情况 下 偏差 是 最 小 的 。 证 明 方 法 是 ， 式 (WU) 取 最 小 时 ， 对 于 A, v 分 别 成 立 ， 
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即将 其 分 别 对 A, v 求 导 应 为 0， 得 


入 十 
fen.) = FO) + FH) -0 ， 

n (53) 
i dvN(v)[x(A,v) - F(A) + F(v)] 20. 
再 做 一 些 代 换 即 可 得 到 式 (GO). 

接 下 来 ， 有 了 式 (bu) 的 形式 后 ， 便 可 结合 式 (bI) 和 (B2) 直接 积分 得 到 某 个 点 (A, v) 的 
五 (7)， 如 此 选取 若干 个 坐标 点 (TACT 程序 选取 了 40 个 点 )， 用 来 三 次 样 条 插值 得 到 所 拟 合 
的 势 ( 样 条 插值 使 拟 合 的 势 会 更 平滑 )。 

(5) 计算 运动 积分 。 借 用 这 个 拟 合 势 ， 以 及 之 前 的 五 , L。， 用 第 三 运动 积分 的 式 (US) 在 
3 个 端点 74 计算 Tg 并 平均 ， 从 而 得 到 3 个 相应 运动 积分 。 其 中 计算 时 需 用 真实 性 条 件 
p2 > 0 进行 判断 ， 如 果 不 符合 ， 就 不 用 端点 而 用 初始 相 空间 点 。 
NS (6) 计算 作用 量 - 角 变 量 。 计 算 的 F(r) 除了 上 一 步 用 来 算 13， 还 用 来 在 式 (GG) 中 以 此 得 
到 当前 点 的 zz， 从 而 计算 作用 量 J 和 和 角 变 量 9 (数值 计算 步骤 具体 见 附录 A) 
2.2.0 Stackel 拟 合 的 其 他 相关 方法 
早 在 1985 年 ，de Zeeuw 和 Lynden-Bell" 就 已 提出 了 把 Stáckel 势 作为 模型 ， 来 局 部 拟 
合 和 全 局 拟 合 星系 引力 势 的 思想 和 方法 。 前 者 在 具体 星系 引力 势 的 平衡 点 附近 展开 引力 势 ， 
然后 通过 对 应 Stáckel 的 椭 球 坐标 形势 来 确定 参数 ， 后 者 通过 以 平均 的 方法 来 近似 Stáckel 
势 的 3 个 分 离 函 数 ， 从 而 计算 运动 积分 。 下 面 简单 介绍 。 

(1) 局 部 拟 合 

在 平衡 位 置 附近 将 V. = V (22,2, 22) 形式 的 引力 势 展开 为 : 


V= Y > > Vietmeryz™ , (54) 


= 其 中 ， 人 另外 ，Stiickel 势 可 写 为 式 (ED) 的 分 离 形式 ， 以 (A, u,v) = 
= (-a, -8, —) 作为 展开 源 点 ， 展 开 分 离 的 函数 选取 在 展开 源 点 附近 的 引力 势 为 0， 即 


Susi 1 二 -1 二 0。 于 是 可 得 : 


32M GO roa", ne) =o n+TB Kv) = mnz+7 . (55) 


anm 


中 用 到 了 (入 + o) 的 展开 式 ， 比 如 : 


Adana {14 y | i yz - y^) | 2y^2 ues] 


b-a y-a (B-a)?  (B—-o)v-o)  (v-ay 


SERAF! 


2 期 


P2 也 同样 地 展开 。 对 比 式 (四 


) 的 形式 ， 式 Ga 


Vooo — 0 , V2oo = —460 , Vaoo = —46i , 


12 
Vi»o = 


BEF, $: Stäckel 引力 势 理论 的 拓展 及 其 应 用 


) 的 展开 系数 应 为 


V220 = 一 一 一 


Veoo = —46» ; (B NH a)? [(G (8 


24[Co(y — B) + mla- y) 4 


a) 


cum O- Ba —3(8 


文献 团 还 讨论 了 这 些 系数 的 约束 条 件 、 存 在 性 
(2) 全 局 拟 合 
文献 回 的 第 4 SORA Sti, Ze HA Shi 

而 且 还 可 以 在 全 局 地 近似 为 Stékel JE 

V — V(z?,y?, 2) 的 全 局 拟 合 方法 。 


式 ， 并 且 给 出 了 椭 球 势 ( 


和 合适 形式 的 选取 。 


圆 星系 有 关 的 引力 势 不 仅 可 以 
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(97) 


局 部 地 ， 


具体 为 镜像 对 称 的 引力 势 
该 方法 要 求 所 拟 合 的 势 在 整体 上 离 Stackel 势 最 近 ， 并 


H24 V 恰好 是 Stäcked 势 时 应 得 到 精确 的 结果 。 下 面 列 出 其 方法 的 主要 步 又 。 
T 对 于 任意 函数 Q = QU, u, v) 定义 下 列 平均 
er Q.A =H J, f QU, u, AMON (v)dudv 
e à demnm 
ie —a ? (58) 
= Q=, ,AN Ma 
N gs MM N (v)dAdpud 
= mun 2 [oa p V)ACA)M QN (v)drdyed 
~ KFP, AA), M(1), N(v) 是 根据 积分 的 约束 选择 的 权重 函数 ， 相 应 的 归 一 化 常数 为 : 
X Ael aja. dr “MGA N=] WG 59 
s [ 40 [mu w= fne (59) 
= 定义 辅助 函数 
XO; i, v) = -(^ — u)y(u — v)(v 和 NV, pv) . (60) 
那么 ， 分 离 函数 可 计算 为 ; 
|. (x-2/3x,  x:(4- (24- n-v) 
io (à — v) 3L 
|. (x-2/3x.0) X- (u— B)DB —v—2) 
_ (x-2/3x,v) Xx: (v-v)Qv-A-np) 
"e sp) 3L 
HH, L- - ig - v)? — X. HERUSCAR p] 中 的 式 (43)， 更 多 的 约束 细节 可 见 该 广 


ChinaXiv& 


162 天 文学 进展 39 卷 


尽管 相 比 于 近年 改进 的 新 方法 (前 述 的 Stáckel 作 浆 法 和 Stáckel 拟 合法 )， 这 些 早期 方 
法 的 适用 性 较 小 ， 但 可 以 说 ， 其 基本 思想 是 相似 的 或 者 被 沿用 了 。 例 如 ， 上 小 节 所 述 的 
Stückel 拟 合法 ， 就 与 这 里 的 拟 合 思路 一 致 : 把 Stackel 势 作为 模型 ， 来 拟 合 具体 星系 的 引 
力 势 。 而 Stáckel E SE1X 5 Stáckel 拟 合法 的 区 别 仅 在 于 : Stáckel 拟 合 法 需要 计算 出 显 式 的 
TAKŽ (7)， 即 局 部 地 基于 某 些 采样 点 把 给 定 的 引力 势 拟 合 为 Stickel 形式 ， 并 计算 分 
离 函数 ， 然 后 对 这 些 采 样 点 上 的 分 离 函数 进行 插值 得 到 最 佳 拟 合 的 Stückel 势 ， 之 后 就 依照 
Stáckel 势 理 论 中 的 各 种 解析 公式 来 计算 运动 积分 及 作用 量 ; 而 Stáckel 作 次 法 是 把 给 定 的 引 
力 势 局 部 地 当 作 Stáckel 势 ， 即 假设 给 定 引 力 势 的 局 部 可 分 离 性 ， 因 此 可 以 采用 一 定 的 技巧 
来 计算 局 部 的 某 种 运动 积分 ， 之 后 就 是 同样 地 依照 Stáckel 势 理论 中 的 各 种 解析 公式 来 计算 
作用 量 。 
2.3 ”基于 环 面 映射 的 方法 

本 节 介 绍 的 内 容 包括 : (1) McGill 和 Binney 提出 、Binney 及 其 合作 者 后 来 完善 的 ， 
通过 正则 变换 构造 环 面 (作用 量 - 角 变 量 ) 来 计算 作用 量 的 方法 ， 即 基于 玩具 模型 的 环 面 映射 
(torus mapping, TM) 思想 及 其 TM PEFR; (2) 基于 此 思想 的 两 个 扩展 : 迭代 环 面 构造 法 
(iterative torus construction, ICTC) ^ — 和 轨道 积分 拟 合生 成 函数 法 (orbit. integration to 
generation function, o2GF) 77, 这 类 方法 的 特点 是 精度 最 高 ， 计 算 花 费 大 ， 属 于 收敛 方 
法 ， 即 可 以 通过 充分 大 的 计算 量 ， 把 作用 量 的 误差 一 直 不 断 地 减 小 下 去 。 

将 环 面 (n-Torus) 视 为 几何 物体 ， 而 将 解析 的 哈密 顿 量 的 环 面 变 形 为 所 求 的 目标 哈密 顿 
量 的 环 面 ， 如 同 微 扰 理 论 ， 通 过 寻求 从 玩具 环 面 坐 标 到 目标 环 面 坐标 的 正则 变换 来 实现 ”。 
即 用 已 知 的 解析 引力 模型 作为 玩具 模型 ， 来 和 逼近 一 般 的 目标 引力 系统 ， 而 这 是 在 环 面 坐标 的 
意义 上 进行 的 ， 并 且 通 过 正则 变换 来 实现 ， 此 即 环 面 映 射 的 关键 思想 。 一 般 所 用 的 玩具 模型 
为 等 时 势 或 谐振 子 势 ， 如 式 (LU) 和 (四 )， 这 些 都 是 Stückel 势 ， 可 解析 求解 。 尽 管 环 面 映 射 
法 不 同 于 前 述 的 通过 近似 Stáckel 的 分 离 函数 F(T) 的 Stáckel 作弊 法 ， 但 其 需要 从 某 种 作 
为 玩具 模型 的 Stickel 势 出 发 ， 本 质 上 相当 于 某 种 简化 版 的 微 扰 法 ， 因 而 仍 可 看 作 是 Stäckel 
势 的 应 用 拓展 。 

本 节 中 ， 若 无 特殊 说 明 ， 用 (J'ev, 0) 代表 玩具 环 面 ( 即 玩具 势 下 粒子 运动 的 环 面 坐 
b). (J, 0) 代表 目标 环 面 ( 即 所 求 的 一 般 引 力 势 下 粒子 运动 的 环 面 坐 标 )。 从 玩具 环 面 映射 
到 目标 环 面 ， 将 此 看 作 一 个 正则 变换 ， 选 取 第 二 类 生成 函数 S = S(T, btey)， 将 这 个 生成 函 
数 对 角 变量 傅 里 叶 展 开 ， 其 一 般 形 式 为 : 

S(9"*, J) = 07 - A(J) -i M S,(J)exp(in- 6°) , (62) 

n0 
其 中 ， 这 里 的 虚 单 位 i 是 为 后 面 求 导 表 达 式 的 简便 而 加 的 。 当 治 着 玩具 环 面 在 0 的 任 一 分 量 
增加 一 个 周期 2r， 相 应 地 ， 目 标 圆 环 上 的 像 点 也 应 完成 一 个 周期 回路 ， 最 终 这 要 求 4 = J. 
于 是 


X 


S(O", J) = gv -J -iX S, (J)exp(in 0*9) . (63) 
n7z0 
通过 作用 量 - 角 变 量 的 实数 性 以 及 哈密 顿 量 的 时 间 反 演 不 变性 ， 可 以 对 生成 函数 的 系数 进行 
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ART: 
Sn =—Sn . (64) 
FEA” 
S(O", J) =O" - J YY S,(J)sinn-0'^", (65) 
nciN 


其 中 ， 向 量 N = {(i, j, k)) 代表 一 系列 精度 的 坐标 (i, j, k) 集 ，(i, j, k) 满足 (Kk > 0) 或 
k=0, j > 0) 3È (k=0, j=0, i 0) 7, xx S(O, J) 求 微分 ， 得 


Os 


toy = = : . pty 
J guy J+ 2 sinan- OY , (66) 
OSs 094,(J) . 
g = — =0 +2 z 0. 
aJ 2 aJ sinn (67) 
这 样 就 通过 生成 函数 联系 起 了 玩具 环 面 和 目标 环 面 。 
> 2.3.1 环 面 映射 程序 


Binney 和 McMillan" $24 Y TM 程序 (Torus Mapper) 站 ， 该 程序 实现 三 个 方面 的 功 
能 ， 有 具体 功能 与 输入 输出 如 下 : 用 户 写 出 目标 引力 势 及 其 一 阶 微分 ， 以 及 给 定 作 用 量 
(Je Jo, J- 1) 对 于 特定 G9， 可 计算 (m, v); 2) 计算 雅克 比 9(z)/6(9); 3) AE (x, J), il 
算 6。 与 传统 的 恒星 轨道 相 比 ， 环 面 的 优点 是 : 一 个 环 面 是 由 100 个 数字 指定 的 ， 即 Sn, 
08,/0J 和 一 些 其 他 参数 ， 而 不 是 由 沿 时 间 序 列 的 数 千 个 相 空 间 坐标 指定 ， 因 此 可 以 在 给 定 
分 状 率 下 显著 缩小 轨道 库 的 大 小 一 。 

环 面 映射 方法 的 流程 见 式 (63)”， 通 过 三 步 正 则 变换 完成 ， 变 换 形式 为 : 


(6, J) 2% (gov, por) 哈 雅 方程 (oy vtoy) OE S) (68) 


已 知 的 真实 势 和 轨道 ， 给 定 目标 作用 量 J 和 系数 S, 的 试探 值 ， 由 玩具 坐标 btey 的 网 格 点 ， 
根据 式 (pa) 可 得 到 相应 的 玩具 作用 量 .Jtey。 在 玩具 势 下 ， 可 直接 通过 解析 的 方法 从 环 面 
坐标 得 到 相 空间 坐 标 (x, utey)。 有 了 这 些 相 空间 坐标 可 计算 相应 采样 点 的 能 量 ， 之 后 用 
Levenberg-Marquardt 算法 欠 代 以 最 小 化 哈密 顿 量 偏差 6 互 ， 得 到 最 终 的 9.。 类 似 地 ， 可 得 
到 905,,/90J， 再 通过 式 (2) 可 得 到 目标 角 变量 o". 

这 样 ， 通 过 真实 的 环 面 坐标 (目标 环 面 坐标 作为 已 知 )， 用 合适 的 系数 S,(J) 和 
905%/0J， 转 化 到 玩具 环 面 坐标 ; 再 通过 玩具 势 下 的 解析 性 就 得 到 相 空 间 坐 标 (xt, vt), 
并 由 点 变换 得 到 (m, wj)。 点 变换 也 是 一 类 正则 变换 ， 在 程序 中 仅 当 /J > 0.05 且 尝 试 在 不 
进行 点 变换 的 情况 下 拟 合 环 面 失败 时 ， 再 使 用 点 变换 一 。 

文献 [EZ] 中 的 附录 图 Al 给 出 了 环 面 映射 方法 的 原理 示意 ， 对 于 一 个 好 的 玩具 模型 ， 更 
高 阶 的 生成 函数 傅 里 时 系数 能 越 来 越 好 地 拟 合 轨道 。 


?TM 程序 的 源码 可 下 载 于 : https://github.com/PaulMcMillan-Astro/Torus。 
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2.3.2 ARABI 

如 上 所 述 ，TM 程序 ”给 出 了 一 般 轴 对 称 系统 的 相 空 间 坐 标 和 环 面 坐 标 转换 的 功能 ， 
主要 是 由 (J, 0) 计算 (x, w)。 对 于 环 面 坐 标的 计算 ， 除 了 按照 上 述 方法 ， 重 复 地 构造 环 面 ， 
直到 得 到 的 相 空 间 坐 标 穿 过 给 定 相 空间 坐标 ，TACT 程序 的 其 中 一 个 方法 对 其 进行 了 改进 ， 
即 迭 代 环 面 构造 法 (TO ^ 7, TIER “Stickel (EME +TM 程序 ”的 迭代 算法 ， 即 用 
Stickel 作 浆 法 计算 试探 的 环 面 坐标 ， 然 后 寻找 环 面 上 与 其 最 近 的 相 空间 点 ， 最 后 补 上 残 差 
FART, 

给 定 一 个 初始 相 空间 点 (z, v)， 先 用 Stückel 作 浆 法 计算 一 个 环 面 坐 标 (Js, 0s): 然后 
用 环 面 映射 得 到 一 个 以 此 为 作用 量 的 环 面 ， 该 环 面 的 频率 为 9， 在 该 环 面 上 找到 一 个 与 给 
定点 最 近 的 相 空 间 点 (ms, vs), HU 


n = |R|?’ (æ — xs)? +(v— vs), (69) 


使 其 最 小 。 之 后 ， 计 算 所 找到 的 这 一 点 的 环 面 坐 标 Joles, vs), Or(as, vs), BS) AJ = 
Jp 一 Js 作为 误差 ， 假 设 这 个 误差 变化 很 慢 ， 于 是 将 其 作为 残 差 补 回去 ， 得 到 一 个 更 好 的 作 
AS J, = Jg + AJ. 
可 以 再 次 构造 一 个 以 Js 为 作用 量 的 环 面 ， 然 后 继续 从 中 找到 一 个 与 给 定 初 始点 最 近 的 
相 空 间 点 ,……， 如 此 往复 迭代 ， 得 到 更 精确 的 作用 量 ， 这 个 过 程 是 收敛 的 。 
2.3.3 轨道 积分 拟 合生 成 函数 法 

这 里 介绍 Sanders 和 Binney™” ”给 出 的 轨道 积分 拟 合 生成 函数 法 。 该 方法 过 程 为 : 轨 
道 时 序 采样 (m, v)(ti) > 计算 玩具 环 面 坐标 (J'07, O°) (t;) — 解 生成 函数 系数 矩阵 。 它 和 
TM 程序 的 区 别 是 ， 该 方法 将 环 面 映射 思想 应 用 于 一 段 运动 过 程 ( 即 一 段 时 间 序 列 )， 计 算出 
在 时 间 t+ 上 沿 轨道 采样 的 一 系列 相 空 间 坐 标的 作用 量 、 角 变量 和 角 频率 ， 即 构造 生成 函数 是 
基于 轨道 积分 ， 而 不 是 TM 所 基于 的 某 个 时 间 点 的 哈密 顿 量 偏差 5 的 最 小 二 乘 。 
具体 步骤 为 ， 给 定 一 个 真实 相 空 间 点 作为 初始 点 ， 首 先 计算 一 个 时 间 Ne T 的 轨道 积分 
( 画 轨道 )， 并 采集 Nam 个 相 空 间 采样 点 ， 其 中 工 是 具有 相同 能 量 的 大 圆 轨 道 的 周期 ， 然 后 
在 每 个 相 空 间 点 ， 解 析 计 算 其 玩具 势 中 的 (ete, .Jt%Y)， 之 后 将 这 些 玩具 环 面 坐标 点 代入 到 
zt (GG) 和 ED 来 求解 整个 方程 组 的 未 知 量 (目标 作用 量 J， 生 成 函数 傅 里 叶 系数 Sn, 及 其 导 
BOS, /OT, UR RA Q 和 常数 O(0)). 

(1) 轨道 积分 与 轨道 分 类 
进行 轨道 积分 。 在 没有 转动 的 情况 下 ， 一 般 三 轴 势 允许 两 类 基本 的 非 共 振 轨 道 : 环 轨道 
和 盒 轨 道 " 一 。 同 一 真实 势 ， 不 同 的 初始 条 件 可 能 导致 不 同类 型 的 轨道 ， 于 是 TACT 为 不 
同类 型 的 轨道 匹配 相应 的 玩具 势 ， 与 给 定 轨道 类 的 环 面 相同 的 几何 结构 。 

(2) 玩具 势 的 拟 合 

对 于 环 轨道 ， 玩 具 势 为 等 时 势 : 


— 


—GM 


Piso = TT oe $ 
b+ VP rr 


(70) 
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对 于 盒 轨道 ， 玩 具 势 为 谐振 子 势 : 


3 
1 2,,2 
Pray p 2 2 Wi Ti : (71) 


其 判定 依据 为 轨道 类 型 ， 即 看 L AL, 的 方向 变化 。 在 用 玩具 引力 势 之 前 ， 需 要 在 轨道 上 
采集 样 点 ， 拟 合 出 相应 玩具 引力 势 的 质量 和 标 长 参数 。 

(3) 生成 函数 

在 每 个 相 空 间 样 点 ti 计算 相应 的 玩具 作用 量 - 角 变 量 (gtey (t), JY) JEFE IERI 
的 式 Ga), HHS J 和 S, KA. Sanders 和 Binney 轨道 积分 拟 合生 成 函数 法 的 实现 技巧 
为 :“ 我 们 不 能 精确 地 解 这 些 方程 ， 因 为 处 理 的 方程 有 无 穷 多 个 未 知 数 。 我 们 可 以 在 每 个 方 
程 的 右边 只 包含 有 限 个 项 ， 则 方程 右边 和 左边 不 能 精确 一 致 ， 因 此 可 以 采取 的 正确 方法 是 : 
对 每 个 方程 的 左右 边 的 差 作 平方 ， 然 后 对 所 有 方程 的 差 平 方 求 和 ( 即 轨道 积分 )， 并 最 小 化 这 
个 总 的 平方 和 ””。 其 总 平方 和 的 表达 式 是 ; 


pu = 2) 2 (ure Joy 人) 一 大 一 2 » ni S4 (J) cos n - se) , (72) 


ncN 


其 中 , k—1,2,3 代表 坐标 维 数 ，N 限制 为 有 限 维 向 量 (有 限 精度 )，|n| < Nma Nmax% 6: 
(ti) 代表 样 点 ， 最 外 层 的 求 和 (关于 指标 i) 是 遍历 所 有 采样 点 ， 此 即 “ 轨 道 积 分 ”的 含义 。 
更 多 关于 样 点 个 数 的 选取 原则 请 详 见 参考 文献 2 中 的 2.2 和 2.3 节 。 

下 面 最 小 化 Ei, BUS EIL, Sm) 的 偏 导 为 0， 


Ji? fe) 
som [4 Yt) — Je — 2 X re Sn(J), cos(n ou) 


ncN 
o= do. -- 233 2m cos (m: el) 。 (03) 
(ro -J2 Y, nS) cos(n 8 ( ») 
ncN 


这 些 方程 可 转化 为 一 个 矩阵 方程 组 


A jroy Å DZ Jtoy = bjtoy 5 (74) 


其 中 ， 
I; cT (t;) 
rtoy = , TJtoy = J,Sn ; 
J p» ( c(t;) elti) . c? (t;) J ( ) 


bioy = S > (J (ts), en (ti) J'Y (t5) , eu (t;) = 2ng cos(n- O(t;)) . 
其 中 ， Cnk (ti) 为 N Fe 3 AY FB 阵 ， TJtoy 和 bto oy 为 (N ae 3) 维 的 列 向 量 ， I; 23 f. pr AB 车 ， 
Ages 为 对 称 和 矩阵 。 最 终 求 解 矩 阵 方程 得 到 目标 作用 量 J 和 生成 函数 傅 里 时 系数 Sno 
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类 似 地 ， 对 于 求 角 变量 8， 则 最 小 化 


-55 (s0 + t; — On (ti )-23. 全 这 san) | (76) 


neEeN 


得 到 角 频 率 Q 以 及 0(0) 和 05,/9J. 
此 外 ， 对 式 (66) 中 的 玩具 作用 量 .JJ'Y 在 玩具 角 变 量 btey 空 间作 平均 ， 得 到 目标 作用 量 


gu. d3gtov 
[e 2 . Btoy 一 iuda NN 
J- ns J 2 SU )sinn | aa (77) 


便 是 一 种 更 简单 的 方法 一 一 平均 生成 函数 法 (average generation function, AVGF) ^ ^ 
2.4 TACT 所 提供 的 方法 的 效果 对 比 

我 们 运行 TACT FU^ ^ 用 几 种 方法 分 别 估算 出 作用 量 ， 如 图 2 所 示 。 我 们 所 使 用 
的 引力 势 为 某 种 轴 对 称 引 力 势 ， 即 McMillan™” 拟 合 的 最 佳 银 河 系 质量 分 布 模型 ， 分 为 薄 
厚 两 个 恒星 盘 、HI 和 Hy 两 个 气体 盘 以 及 核 球 和 暗物质 晕 ; 各 个 星系 成 分 的 模型 输入 参数 
(适合 TACT 程序 的 格式 )， 已 在 作者 的 github 主页 加 给 出 。 我 们 采用 的 恒星 轨道 为 典型 的 
银河 系 薄 盘 轨 道 ， 初 始 条 件 为 (zx, v) = (8.178, 0.1, 0.225; 20.13，187.01，14.95)， 单 位 为 
kpe 和 kmy/s。 图 中 的 几 种 方法 分 别 是 : Fudge v1 和 Fudge v2 为 Stáckel Fudge Ek, H 
rH, Fudge vl 用 式 (HN) 估算 焦距 ，Fudge v2 基于 壳 轨 道 (shell orbits) 来 估算 焦距 ，Fit 对 
应 Stückel 势 拟 合法 ，ItTC 对 应 环 面 迭代 法 ，O2GF 对 应 轨道 积分 拟 合生 成 函数 法 ，AvGF 
为 平均 生成 函数 法 CAA 和 SAA 分 别 为 柱 绝热 近似 和 球 绝热 近似 方法 ， 下 一 节 作 简介 。 轨 
道中 的 每 个 点 都 独立 计算 一 次 作用 量 ， 因 而 作用 量 的 数值 随时 间 有 所 波动 。 


= CAA —SAA 

= Fudge vl == Fudge v2 
= Fit =" ItTe 

== AvGF == O2GF 


. sol = 
J,/kpe-km's J,/kpe-km-s™ 


图 2 ”作用 量 数值 计算 的 效果 对 比 


Sanders 和 Binney ”对 几 种 方法 的 在 轨道 类 型 、 精 度 和 时 间 等 方面 的 性 能 做 了 更 详细 
的 对 比 ( 见 其 中 的 图 2 到 图 6)。 总 体 上 看 ， 收 敛 方 法 比 非 收敛 方法 有 更 高 的 精度 ， 但 耗 时 也 


高 1 ~ 3 个 量 级 。 其 中 ， 轨 道 积分 拟 合生 成 函数 法 精度 最 高 ， 耗 时 较 长 ， Stückel 作 浆 法 速 


https: //github.com/PaulMcMillan-Astro/GalPot /blob/master/pot/PJM17_best.Tpot. 
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度 很 快 ， 是 大 量 计 算 的 最 经 济 方法 ， 但 精度 不 如 收敛 性 方法 高 ， 特 别 是 对 于 盒 轨 道 的 误差 
大 。 另 外 ， 对 于 轴 对 称 势 经 过 环 面 迭代 后 (BU ICTC 法 ， 目 前 TM 程序 只 提供 了 轴 对 称 的 环 
面 迭代 )，Stéckel 作 次 法 的 精度 会 提高 ， 而 耗 时 也 变 长 。 此 外 ，TACT 程序 包 中 ， 轨 道 积分 
拟 合生 成 函数 法 和 Stáckel 作 浆 法 还 能 用 于 三 轴 势 。 
2.5 ”其 他 一 些 求 运 动 积分 的 方法 或 程序 
前 面 两 节 我 们 主要 介绍 了 Binney 及 其 合作 者 们 近年 来 的 工作 。 接 下 来 我 们 介绍 其 他 一 
些 基于 Stáckel 势 的 估算 运动 积分 或 作用 量 的 方法 。 
2.5.1 “近似 的 运动 积分 ”解析 式 
Bienaymé 给 出 了 一 种 通过 椭 球 坐标 轴 上 的 代表 点 来 计算 所 有 位 置 处 的 引力 势 的 方法 ， 
因此 一 般 引 力 势 可 以 拟 合 为 Stackel 形式 。 然 后 ， 就 可 以 根据 Stackel 势 的 运动 积分 解析 表 
达 式 ， 给 出 这 个 一 般 引 力 势 下 的 “近似 的 运动 积分 ”解析 式 (又 称 为 “ 准 运动 积分 ”，quasi 
integral of motion" ， 特 别 重要 的 是 “quasi- 厂 ”的 表达 式 )。 该 方法 适应 接近 Stickel 形式 
的 势 。 下 面 介 绍 这 种 方法 。 
r 由 式 (ZI), 将 PRA, WSN: 


p- 4.8090 * o90 Y 80 t v) _ nulu + a)l + Bn t v 
(^ — p)(A- v) (u — v)(u — A) (78) 
SOV tay t+ Bv+y 
(v — X)(v ~ n) l 


3k o f£ (—o, —8, —y) 连续 ， 应 有 Ba, —B, —y) = 0。 让 势 两 边 取 (A, —8, -7) 位 置 
的 值 ， 则 有 
N $(A, —8,— 
ae ca) = - SRM), (79) 


qu 同 理 有 


= _(—a,p, =) xui es Em —p,v) 
= id (neg ^ "0 A(v +7) 0) 
将 上 面 的 分 离 函 数 代 回 到 式 (ZB), AAS 
(A B) t 9) (K+ o)(u a | 
THES Dew m e zo y* , 1, —-7)4 m 
(v+ OUE 5 PEN 
(v — A)(v — u) s 


ROME RAED AS "AERA ER" (Uu (A, 一 6, —y) 中 后 两 个 坐标 固定 ) 上 的 势 值 来 表示 任 
意 点 的 引力 势 。 
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采取 类 似 的 手法 ， 这 个 引力 势 下 的 运动 积分 可 以 解析 性 地 表达 出 来 。 运 动 积 分 J, K 的 
表达 式 为 
J =1/2L? — 1/2(8 + y)v? — 1/2(y + oa)w — 1/2(a + 8)v2 


ur- (922 +) Sy 


K =1/2B-yv2 + 1/2yov? + 1/2o8v2 — 1/2aL2 — 1/282 — 1/2412— 
C(A 
" o 2 m 
相关 推导 见 参考 文献 回 中 的 式 (9) 和 (14)。 将 式 (Izd) 和 (BO) 代入， 可 计算 得 到 与 式 Cg) 同 
样 方式 的 运动 积分 J, K 的 表达 式 。 
类 似 地 ， 由 式 (15), D, Is (第 一 个 运动 积分 是 能 量 E 


È E) 的 表达 式 可 写 为 
I = H(A, p, v) 1/2(a Bv? | V2 


= Ij = € p, v) +1/2(y — 8)v? 4 Ar + 1/202 
即 参考 文献 BA 中 的 式 (9) 和 (10). oci [ba] 推导 了 上 式 中 的 椭 球 坐标 的 函数 项 ， 得 到 
Ge (34) 
E= egy. v) , 
KH, BO, p,v) 即 为 式 (I8) 中 的 形式 ， 其 简化 形式 见 文 献 Bal] 中 的 式 (16) 和 (17). 
该 方法 计算 “每 一 点 的 ”运动 积分 ， 需 要 当前 点 的 坐标 、 速 度 和 角 动 量 等 轨道 信息 。 文 
献 [ba] 在 后 文 用 高 度 精确 的 数值 方法 计算 出 轨道 ， 并 研究 了 这 种 方法 所 计算 的 运动 积分 的 
误差 ， 最 后 做 了 Jeans 方程 的 检验 。 
2.5.0 ”绝热 近似 
对 于 轴 对 称 势 ， 文 献 [22] 还 介绍 了 用 近似 公式 来 估算 作用 量 的 方法 。 通 过 假设 在 柱 坐 
标 系 或 球 坐 标 系 中 引力 势 是 可 分 离 的 ， 可 求 出 这 些 作用 量 。 
(1) 柱 绝热 近似 ”。 这 种 观点 认为 ， 得 轨道 恒星 的 重 向 (HI 2 方向) 频率 OO, 比 其 径 向 ; 
率 Q. 大 得 多 ， 因 此 当 恒 星 径 向 振荡 时 ， 主 导 其 垂 向 振荡 的 势 可 被 认为 变化 缓慢 ， 其 结果 误 
是 是 绝热 不 变 的 ”。 用 这 些 假设 推导 出 作用 量 的 近似 公式 


g jem 
太志 一 dzv, 
T Jo 


1 ffe | 
jte ri dRvg 


T Rperi 
以 及 计算 积分 上 下 限时 需要 用 到 的 垂直 能 量 ( 即 垂直 动能 加 垂直 势能 ， 其 中 前 者 是 由 于 速度 
的 平方 按 正 交 坐 标 方向 分 为 三 个 分 量 平方 的 相 加 ; 后 者 是 保守 力 按 照 某 方向 的 分 解 所 得 的 


2 2 
I2 vin 
(83) 


y= 


SS 
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势能 ) 和 有 效 势 : 
E jud + &(R, z) — 6(R,0) 
2 , (86) 


Js 
Gr = ó(R,0) 十 Dd + Bs(Js, R) — BE,(J;, Re) 


RP, Re 是 导 心 半径 (可 参见 有 关 “ 本 轮 近 似 ” 的 文献 )。 
(2) 球 绝热 近似 假设 治 长 球 坐 标的 球面 存在 绝热 不 变 振荡 。 作 用 量 的 近似 公式 为 : 


2 Vy 9 Vy 
J — = | dvp, = 7/ dv /2P2(A, v) / E, E $, (v) 


1 p% "EY ; (87) 
d =f dApA = =f d\y/2P2(A, e) / Ex — (3) 
其 中 ， ; 
$,(v) = A uc B(A, v) — BA, c?) 
S, (A) = 90,08) + + ER = VITA, Jg, J) (88) 


具体 计算 可 见 参考 文献 [24]. 
2.5.3 ”一些 相关 程序 

除了 Sanders 和 Binney” 的 TAT 程序 包 以 及 其 环 面 迭代 功能 所 依赖 的 TM 程序 包 ， 
星系 动力 学 还 有 一 些 程序 可 供 使 用 ， 这 些 程序 现在 都 可 以 公开 获取 。 主 要 程序 如 下 所 述 。 

(1) Gadget, Gadget 以 及 其 后 续 发 展 Arepo8 程 序 ， 是 “引力 N 体 + 流 体力 学 ”模拟 
程序 。 在 星系 动力 学 领域 ， 我 们 可 以 用 它 来 做 无 碰撞 N 体 模拟 (特别 是 当 我 们 考察 暗物质 粒 
子 构成 的 暗物质 晕 的 动力 学 时 )。 

(2) Galpy，Galpy (a ud library for galactic dynamics)™ 是 Bovy " 提供 的 Python 
程序 包 ， 能 用 来 进行 数值 轨道 积分 、 作 用 量 - 角 变 量 、 分 布 函数 等 相关 计算 。 

(3) AGAMA, AGAMA (action-based galaxy modelling architecture)E! 是 Vasiliev" #2 
供 的 星系 动力 学 程序 ， 能 提供 任意 解析 密度 剖面 或 N 体 模 型 的 引力 势 的 计算 、 轨 道 积 分 和 
分 析 、 位 置 -速度 与 作用 -角度 变量 之 间 的 转换 、 分 布 函数 的 计算 以 及 迭代 构造 自治 多 组 分 星 
系 模型 等 功能 


M 


° https: //wwwmpa.mpa-garching.mpg.de/galform/gadget /index.shtml 
?https://arepo-code.org 
?https://github.com/jobovy/galpy 


?https://github.com/GalacticDynamics-Oxford/ Agama 
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2.5.4 深度 学 习 在 作用 量 计 算 中 的 应 用 

在 这 一 小 节 中 ， 我 们 介绍 最 新 的 一 个 研究 动态 。Ibata 等 人 ”将 人 工 智 能 (AT) 中 的 深 
度 学 习 方法 应 用 到 作用 量 计算 上 。 该 文 用 AI 算法 来 实现 在 静态 引力 势 下 从 粒子 在 相 空间 
的 “位 置 -速度 ”到 “作用 量 - 角 变 量 ” 的 映射 ， 并 且 得 到 加 速度 场 。 该 算法 的 核心 原理 与 前 
述 “ 生 成 函数 法 ”类 似 ， 先 从 数据 点 中 拟 合 玩 具 势 ， 然 后 通过 从 玩具 环 面 (9, J) 到 目标 环 
面 (0, J’) 的 正则 变换 表达 式 ， 和 迭代 出 最 合适 的 生成 函数 G; 之 后 是 正则 点 变换 的 生成 函数 
P， 最 终 得 到 所 求 的 目标 环 面 (9”, J”)。 其 中 两 个 正则 变换 所 用 网 络 的 形式 为 : 


G = netg(cos 0, sin, J’) , P —netp(cosO, sin0) . (89) 


该 算法 的 程序 为 ACTIONFINDER， 基 于 Python 的 Pytorch 所 写 。 据 文献 BQ] 介绍 ， 
经 测试 ACTIONFINDER 得 到 的 作用 量 精 度 要 比 Stickel 作 次 法 要 高 ， 而 且 它 还 有 如 下 优 
点 : 相 比 于 早期 的 Torus Mapper 对 初 值 要 求 更 不 敏感 ， 不 需要 向 它 提供 系统 内 轨道 的 公平 
样本 ， 几 乎 可 履 盖 任何 感 兴趣 区 域 的 样本 ; 不 需要 预先 知道 所 研究 系统 的 哈密 顿 量 或 引力 


- 势 ， 此 外 ， 该 算法 所 用 的 深层 神经 网 络 ， 相 比 前 述 “ 生 成 函数 法 ”线性 分 解 成 的 傅 里 叶 系 
> 数 ， 能 够 生成 更 加 灵活 的 非 线性 函数 ， 并 且 这 种 灵活 性 可 能 使 其 更 容易 拟 合 更 一 般 的 动力 系 
统 。 看 来 ， 将 来 经 过 一 定 的 发 展 完善 后 ， 深 度 学 习 或 其 他 AI 技术 用 来 计算 作用 量 ， 适 用 性 
会 比 传统 方法 更 广泛 、 性 能 更 佳 ， 应 该 有 一 定 的 发 展 前 景 ， 未 来 或 成 为 分 析 更 复杂 体系 的 另 
一 种 尝试 。 
3 ”应 用 


在 星系 研究 中 ， 作 用 量 的 一 个 重要 应 用 在 分 布 函数 方面 ， 即 使 用 以 作用 量 作为 自 变量 的 
分 布 函数 ， 为 星系 建 模 ”。 一 般 来 说 ， 我 们 难以 直接 获得 分 布 函数 的 形式 。 在 建 模 以 及 数 
据 拟 合 的 实际 工作 中 ， 人 们 通常 是 为 不 同 的 星系 组 分 (如 恒星 盘 、 核 球 、 暗 物质 尝 等 ) 采用 
某 种 先 验 的 分 布 函数 形式 ， 然 后 根据 数据 拟 合 来 确定 其 中 的 具体 参数 的 取 值 。 

KY 10 年 前 ， 基 于 分 布 函数 的 建 模 方法 实际 上 很 难 应 用 ( 见 Read ”和 Binney” 的 综 
= 述 )， 星 系 的 建 模 工作 往往 是 基于 Jeans 方程 ， 即 对 分 布 函数 所 满足 的 无 磁 撞 玻 尔 兹 曼 方 程 
[T BOR (BU ATI “RE ATI” ). [BoERRARTTSSS ( 即 OF /Ot = 0)， 从 无 碰撞 玻 尔 兹 曼 方 程 出 发 ， 
分 别 乘 上 速度 分 量 并 且 积分 掉 速 度 空间 ( 即 取 和 矩 )， 即 得 到 Jeans 方程 。 对 于 实际 的 银河 系 建 
模 工 作 ， 研 究 者 还 会 假设 轴 对 称 ( 即 在 柱 极 坐标 下 9f/09 = 0, Of /Ovg — 0)， 那 么 Jeans F 
程 的 表达 式 为 ( 柱 极 坐标 ): 


OR Oz R OR 

1 0(R^vvmgv;) | O(vv;v;) | ; (90) 
RP ØR bus ona dS 

1 O(RvoR.)  O(vo , Ob . i 

HR OH ' ws "ga $e 90M 


7V1 


a 
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其 中 , v 为 示 踪 者 恒星 的 数 密度 ， 为 分 布 函数 的 0 阶 矩 ;平均 速度 元 为 1 阶 和 矩 ， 速 度 弥散 张 


o1 
^v 


为 2 阶 和 矩 。 在 观测 上 ， 速 度 弥散 可 由 某 一 个 位 置 处 诸 恒星 的 速度 方差 来 计算 。 

正如 Read MAET, Jeans 方程 和 矩 方法 的 优点 是 : 与 其 他 方法 相 比 非常 快 ， 可 以 探 
索 较 大 的 参数 空间 ;并 且 不 需要 假设 分 布 函数 的 形式 ， 因 为 只 是 限制 了 它 的 矩 。 缺 点 是 : 必 
须 对 数据 进行 整理 以 计算 速度 的 2 RE (但 实际 上 恒星 的 速度 数据 往往 不 是 高 斯 分 布 ， 因 此 
不 能 完全 为 二 阶 矩 所 刻画 ); 分 布 函数 的 形式 不 能 使 用 (因此 未 能 使 用 恒星 完全 的 分 布 信息 ); 
Jeans 方程 组 一 般 情况 下 不 能 闭合 ( 即 不 能 求解 )， 在 某 些 情况 下 可 能 得 到 一 个 实际 分 布 函 数 
不 存在 的 解 。 

如 果 这 时 候 我 们 进一步 假设 该 星系 中 的 恒星 运动 只 遵守 两 个 运动 积分 ( 即 文献 中 常 表述 
的 f(E,L.). WA or = 0, M@Hor, =0. WN, z 方向 的 Jeans 方程 ， 就 仅 与 各 量 的 z 
方向 的 梯度 有 关 ， 与 RR 方向 的 运动 解除 看 合 了 ， 即 : 

2 
ee ^ 

在 给 定 v, $ 之 后 ， 我 们 可 以 得 到 op (或 os) 的 表达 式 以 及 吧 的 表达 式 ， 也 就 是 说 Jeans 7j 
程 闭合 了 ”。 

在 更 早 之 前 ， 就 连 上 述 2 个 运动 积分 情形 下 的 Jeans 方程 都 难以 得 到 银河 系 观测 数据 
(恒星 的 密度 分 布 数据 以 及 运动 学 数据 ) 的 支持 ， 研 究 者 把 引力 势 的 表达 式 ( 即 泊 松 方程 ) 也 
“一 维 化 ”了 ， 简 化 为 : 


Cij 


J Eo- m-e), (91) 


9?ó 
— 
zt (E2) 和 (E3) 构成 了 所 谓 的 “一 维 近似 ”加 。 由 于 —09/8z WE K (引力 的 z 方 向 分 
量 )， 早 期 文献 中 又 把 这 种 一 维 近似 称 作 “K, force” 方 法 ”。 如 果 用 更 方便 的 观测 量 恒星 
面 密度 (X) 来 代替 密度 p， 那 么 泊 松 方程 可 以 表示 为 K, = 2nGd. 

在 下 面 的 几 个 小 节 中 ， 本 文 介 绍 近年 来 有 代表 性 的 建 模 工作 。 这 几 个 工作 都 是 基于 3 
个 运动 积分 (或 作用 量 ) 的 分 布 函数 来 为 银河 系 建 模 。 这 些 工 作 中 ， 对 于 观测 到 的 恒星 盘 成 
分 ， 一 般 是 使 用 基于 作用 量 的 准 等 温 球 分 布 函数 模型 ， 或 者 基于 运动 积分 的 准 等 温 球 分 布 函 
数 。 对 于 暗物质 党， 这 些 工 作 使 用 质量 分 布 模型 (如 NFW 模型 ) 而 不 是 分 布 函数 模型 。 后 
K, Pim ATK Binney 和 Pi 组 开始 尝试 使 用 分 布 函数 模型 来 直接 为 暗物质 晕 建 模 并 拟 
合 银河 系数 据 ， 这 项 暗物质 晕 建 模 工 作 由 于 数据 有 限 以 及 软件 上 的 问题 ， 目 前 的 效果 不 太 理 
想 宇 。 在 本 节 的 末尾 ， 我 们 还 介绍 前 人 提出 的 针对 恒星 盘 的 另 一 种 分 布 函数 模型 。 

前 人 的 这 些 工作 ， 主 要 是 基于 RAVE 和 SDSS 的 银河 系 巡 天 数据 。 现 在 ， 在 Gaia 的 天 
测 和 测 光 数据 以 及 大 规模 光谱 巡天 (例如 中 国 的 LAMOST) 的 基础 上 ， 我 们 预计 这 些 基 于 分 
布 函数 建 模 的 方法 将 大 放 异 彩 。 


AnGp = (93) 
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3.1 BF f(Jn, Le, J.) 的 银河 系 动力 学 建 模 
3.1.1 Bovy 和 Rix 的 模型 

Bovy 和 Rix” 用 SEGUE 的 G 矮星 样本 (范围 为 5 < Rac < 12kpc, 0.3 < |z| < 3kpc, 
HEP, Rae 为 到 银河 系 中 心 的 径 向 距离 ，|z| 为 到 银河 系 中 心 的 垂 向 距离 ) 进行 建 模 ， 来 
推断 恒星 和 暗物质 质量 分 布 。 将 数据 按 金 属 丰 度 分 为 43 个 单 丰 度 亚 群 (mono-abundance 
population, MAP)， 对 每 个 MAP 的 相 空间 数据 独立 地 设置 一 个 参数 化 的 银河 系 引 力 势 和 一 
个 基于 作用 量 的 准 等 温 的 恒星 分 布 函数 ， 并 且 假 设 每 个 MAP 的 数据 在 不 同 的 半径 约束 垂 向 
力 为 尺 .。。 用 这 些 结果 ， 研 究 了 在 4.5kpc < Rac < 9kpc 范围 ， 银 河 系 中 平面 1.1 kpc 内 的 
总 面 密度 D1 (R) = ER dzp(R, z) 和 重 向 力 Ko VAT ERS S PRAE PLA 

(1) 恒星 的 分 布 函数 模型 

准 等 温 分 布 模型 (quasi-distribution function, gDF)” 是 一 种 最 常见 的 恒星 分 布 模型 。 
将 示 踪 星 的 分 布 函数 建 模 为 一 个 或 多 个 等 温 分 量 ， 根据 定义 ， 等 温 分 量具 有 恒定 的 速度 色 
B tz 无 关 ， 从 而 速度 分 布 必须 是 高 斯 分 布 ”。Bovy 和 Rix™ 文中 的 每 个 MAP 都 是 由 
准 等 温 分 布 拟 合 的 恒星 种 群 ， 分 布 函数 的 自由 参数 ， 与 示 踪 剂 群体 的 径 向 密度 分 布 ， 径 向 和 
垂 向 速度 色散 ， 以 及 它们 如 何 随 RR 变化 有 关 。Bovy 和 Rix 宇 所 用 的 基于 作用 量 的 恒星 分 布 
模型 来 自 文 献 6]: 


V vJ, 
qDF (a, v) = forn(Jr, Lz) x 2no2 x exp (xx) ; (94) 
Rn( Re) KJR 
Rr, forl(Jr, Lz) = mox x [1 + tanh(L;/ Lo)] x exp “Re } k, Q, v 是 本 轮 


频率 、 贺 周 频 率 和 垂 向 频率 ，RR。 是 具有 和 角 动 量 L 的 圆 形 轨道 的 半径 。 并 且 假 设 示 踪 者 密度 
径 向 轮廓 、 径 向 和 垂 向 速度 弥散 是 指数 形式 的 : 


n(R.) x exp(— R./hg) 
or(R.) = oro exp(—(R.— Ro)/ho,) -> (95) 
c, (Re) = 02,0 exp(—(R. — Ro) /he.) 


分 布 函数 中 的 作用 量 ， 采 用 文献 [E] 的 Stäckel 作 浆 法 来 计算 。 
(2) 银河 系 引力 势 模型 
为 计算 轨道 作用 量 ，Bovy 和 Rix 四 提出 一 个 完整 的 三 维 模型 来 计算 银河 系 的 势 ， 一 个 
先 验 有 许多 自由 参数 的 “ 核 球 + 恒 星 盘 + 气 体 盘 十 暗物质 晕 ” 模 型。 
指数 盘 密 度 为 双 宕 律 盘 模型 ， 


Paisk(R, z) = pa exp( — R/ Ra — |Z|/z4) ， (96) 


引力 势 通过 该 密度 的 泊 松 方程 的 解析 解 得 到 。 


pbulge(7) = d E m ; 
(r/a)*(1 + r/a(8 — a) 


(97) 
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其 中 ，aw = 1, 6 = 4 时 ， 称 为 Hernquist 模型 ，a —2, 8 =4 时 ， 称 为 Jaffe 模型 。 暗 物质 晕 


AE Gp. 


Pom(R, z) e 1/r* , (98) 
它们 的 共同 引力 势 为 总 引力 势 。 
在 Bovy 和 Rix” 对 单个 MAP 的 qDF 拟 合 中 ， 除 了 恒星 盘 的 标 长 Ry MEME Ro 
处 的 径 向 力 的 贡献 以 外 ， 其 他 基本 参数 (TEBE zr LAR Ro 附近 圆周 速度 ve(Ro) 和 
旋转 曲线 的 “平坦 度 ”dlgvc(Ro)/dlg(R)) 都 被 固定 。 
(3) 拟 合 手续 
Bovy 和 Rix ^ 所 使 用 的 拟 合 方法 为 最 大 似 然 法 ， 其 拟 合 步骤 如 下 ， 流 程 示意 图 见 图 3。 


在 参数 p,, p, 所 构成 空间 的 每 一 点 ; 


示 踪 者 恒星 观测 数据 
(1, b, D, V, rng [Fe/HI]) 


起 于 观测 量 的 恒星 密度 = 
AU, b, D, v, r, gr, [Fe/H]lp,, py) | | 银河系 信息 


EHEJ T 
I 
单个 数据 点 的 恒星 概率 
qDF(J[x, v]lpy. Por) 
最 大 似 然 法 


BRAKA 


3 Bovy 和 Rix MA 合流 程 的 示意 图 


观测 量 下 的 分 布 为 : 
= A(L b, D, v, r, g — r, [Fe/H]) 2p(r, g — r, [Fe/H] ||R, z, 9) x qDF (a, v)x (99) 
全 |J(R, z, p; 1,6, D)| x S(plate,r,g—r) , 


其 中 ，D 为 观测 量 距离 ，! 为 银 经 ，&b 为 银 纬 ，g 一 7 ABH, [Fe/H] 为 观测 量 金属 丰 度 ，v 
为 速度 ， 由 观测 得 到 ，(R, 2,6) 为 银 心 柱 坐 标 系 ， 坐标 变换 的 雅克 比 ，5 为 选择 函数 。 
归 一 化 后 ， 每 个 数据 点 的 似 然 为 : 


In.Z; pp — InqDF(J(zi,vi)|pa, ppr)— 


(100) 
m f dldbdDdvdrd(g — r)d[Fe/H]A(l, b, D, v, r, g — r, |[Fe/H]|ps, por) . 


加 上 参数 条 件 是 表示 设 参数 po, por 取 某 一 组 值 的 条 件 ， 最 大 似 然 就 是 在 参数 取 值 空间 中 
找到 似 然 最 大 时 的 参数 值 。 作 用 量 的 计算 会 用 到 总 引力 势 8(R,z)， 这 些 引 力 势 的 参数 实际 
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是 待定 的 基准 势 ; 正 是 在 给 定 分 布 函数 模型 和 引力 势 模型 后 ， 用 最 大 似 然 法 依靠 数据 从 参数 
空间 找到 最 合适 的 参数 后 ， 才 能 得 到 最 合适 的 模型 。 因 此 ， 引 力 势 的 参数 ， 正 是 通过 数据 对 
作用 量 的 依赖 ， 输 入 到 观测 到 的 密度 当中 。 

数据 点 是 独立 的 ， 考 虑 异常 数据 点 (outliers) 后 ， 总 的 似 然 为 : 


In.£ — Y ^ In[(1 — pos) In Zip + pou In Gout] - (101) 


最 后 根据 最 大 似 然 法 估计 ， 使 这 个 总 似 然 最 大 ， 即 找到 所 求 参数 po, por 的 最 合适 值 。 
(4) 结果 
利用 这 些 模 型 的 所 有 参数 的 最 佳 结 果 ，Bovy 和 Rix 得 到 了 星系 盘 的 径 向 标 长 Ru = 

(2.15+0.14) kpc， 与 通过 光度 法 推断 的 恒星 质量 的 空间 分 布 一 致 。 

来 自 不 同 的 MAP 数据 ， 以 不 同 的 半径 约束 a(R): 因而 Bovy 和 Rix 用 这 43 个 
MAP 数据 的 平均 半径 和 Y. 的 平均 值 的 关系 ， 研 究 了 2 1 随 半 径 的 变化 ， 这 些 MAP 的 结 
果 也 一 致 ”。 类 似 地 ， 还 得 到 了 1.1kpc 的 垂 向 力 Kaia (R) AASR El 距离 标 度 的 后 
验 分 布 函数 所 导出 的 最 佳 半 径 ， 将 这 些 分 别 拟 合 为 : 

| 51.1(R) = 70Mopc ^? exp(—(R — Ro)/(2.5 kpc)) 
K,ia(R) = 68Mopc ? exp(—(R — Ro)/(2.7 kpc)) 


(102) 


Bovy 和 Rix ^" 在 计算 恒星 盘 对 圆周 速度 的 贡献 、 暗 物质 晕 的 寡 律 指数 等 结果 方面 ， 对 
完整 Jeans 方程 和 垂 向 力 两 种 方法 进行 了 对 比 ( 见 文 献 [Ha] 中 的 图 18, 19, 20); 还 展示 了 恒 
星 盘 和 暗物质 尝 分 别 对 旋转 速度 的 局 部 贡献 ( 见 文 献 (a) 中 的 图 21); 得 到 的 对 暗物质 军 的 
径 向 轮廓 的 限制 为 : 在 95% 置信 度 下 ，ppxw(r; URO) x 1/r^, a < 1.53. 

3.1.2 Pif 等 人 的 模型 

与 Bovy 和 Rix 提出 的 模型 类 似 ， 同 样 用 数据 限制 银河 系 的 参数 ，P 进 等 人 用 距 银 
河 系 平面 约 1.5 kpc 内 约 200000 颗 巨 星 的 运动 学 ， 来 测量 太阳 附近 质量 密度 的 垂 向 分 布 ， 
以 及 局 部 暗物质 嘿 贡 献 。 其 所 用 的 模型 的 组 分 更 多 ， 由 气体 盘 、 薄 盘 、 厚 盘 、 核 球 和 暗物质 
深 五 个 部 分 组 成 ， 下 面 介绍 其 质量 模型 和 分 布 函数 。 


(1) 质量 模型 
盘 的 密度 模型 为 : 
Xo R Z Rhoe 
R = 103 
pi 2) = exp| (zie R ) , ( ) 


其 中 ， 气 体 盘 非 零 参 数 Rio 是 为 了 在 圆 盘 中 创建 中 心 腔 ， 而 其 他 两 个 盘 则 设置 为 零 。 
核 球 和 暗物质 晕 的 密度 模型 为 : 


p2 (R, z) = 


and may Pl (mro/rent)’] ， (104) 


R ? 之 a ] LM = SL, MA ZA 一 AAA ae 
Hob, m(R,z) = | (=) + (=) ，po 设置 密度 标 度 ，7o 为 半径 标 度 ，g 为 等 密度 面 的 
0 
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ae 


HLL, cut 为 截断 ， 指 数 y, B 控制 径 向 密度 轮廓 的 内 外 径 向 斜率 ; TRR, AT 
别 采用 不 同 的 轴 比 和 截断 等 ”。 公 式 中 的 其 他 参数 的 意义 与 上 一 节 的 相同 。 

因此 Pim 等 人 的 质量 模型 具有 8 个 参数 ， 注 盘 和 厚 得 分 别 具 有 3 个 自由 参数 
(Xo, Ra, za). MH EROR 2 个 自由 参数 po,ro， 其 中 薄 盘 和 厚 盘 的 径 向 标 长 取 相 
IR]. SCHR (BO) 中 的 表 1 给 出 了 银河 系 模型 的 固定 参数 取 值 。 

(2) 分 布 函数 模型 

分 布 函数 分 为 盘 和 尝 两 大 部 分 ， 


f(D, Jz; Lz) = Pais Ss J 55 Lz) + PUN Jr, J xs Lz) , (105) 


RP, fas UERITAS AB, fhao SERDARA. Fao AES EIU CZ Ls 
盘 的 分 布 函数 为 ; 


faisk (i, Ja Lz) = finr Jay Lz) + Fynick fthick (Jr; Jay Lz) , (106) 


RP, finn 为 注 盘 的 分 布 函 数 ，ftnia 为 厚 盘 的 分 布 函数 ，Finick 是 厚 盘 与 沙盘 的 质量 之 比 。 
盘 的 分 布 函数 也 是 基于 作用 量 的 ， 并 且 通 过 Stückel 作弊 法 ” 估算。 类似 于 前 人 的 准 等 
温 模型 ， 盘 分 布 函数 也 是 由 准 等 温 组 分 构成 ， 


fo (J Ly) = zu + tanh(L/ Lo)| exp(—&J../o2) 


fo, (Ji, Lz) = ane gg vJ; /02) 


Xp. v 仍 代表 垂 向 频率 而 不 是 椭 球 坐标 。 但 是 Pi 等 人 ”的 模型 中 ， 注 盘 还 考虑 了 星 族 
年 龄 ， 其 速度 弥散 除了 工 ， 还 依赖 于 年 龄 r， 


D(L) = Xoexp(—R./ R4) 


TTT f 
Ara cpu PSA (108) 
B 
THT 
o,(Lz,T) = 04 (=>) exp(Ro — Re)/Ro, 


3t — HW AB Ta P EER BERE T] 8 2C P E, REIER EREA to TEE R 
TU: 


P dre™/* fo, (ds La) fo, (Jz, Lz) 
(e7m/to = 1)to 
Pit 4& AP" vr fe 4) ACI EN fro UB LED er E Ee 
ALIS, MATUAEEEREAe ARI IEEE. ERNARI Posti 等 全 REA HUS 


an(J Jz, Dz) — (109) 


fhalo (J, Jz, Lz) = (A) ; (110) 
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12, TNT: 
Rm, h(I) = J, + o^ 十 万 Ies h(J) 是 J 的 一 次 齐 次 函数 。 


Pim 等 人 采用 最 小 二 乘法 来 建立 这 些 模 型 ， 先 给 出 了 最 佳 拟 合 的 待定 参数 ， 然 后 给 
出 了 各 组 分 的 密度 轮廓 ， 他 们 还 分 别 给 出 了 重子 物质 与 暗物质 曼 对 圆周 速度 和 中 盘面 总 面 
密度 随 径 向 的 局 部 贡献 变化 。 
3.2 EF f(E, Lz, i) 的 银河 系 动力 学 建 模 
前 面 几 个 准 等 温 模型 及 其 变形 都 是 基于 作用 量 的 分 布 函数 模型 。 除 了 基于 作用 量 的 ， 以 
前 也 有 直接 基于 E, La I5 三 个 运动 积分 的 分 布 函数 模型 。 关 于 第 三 运动 积分 ， 大 体 上 也 是 
以 Stackel 引力 势 中 的 运动 积分 的 形式 为 蓝本 ， 作 近似 假设 。 下 面 简 介 其 中 的 一 些 模型 。 
3.2.1 Bienaymé 等 人 的 模型 

Bienaymé 等 人 ”也 用 RAVE 的 红 团 复 星 数据 研究 了 径 向 力 和 暗物质 晕 的 局 部 贡献 。 其 
假设 银河 系 引 力 势 为 某 种 Stackel A: 


111 
Y À —y ? ( ) 


Heh, h(r) = GM /T q H Kuzmin-Kutuzov $5" "36 (多 了 参数 g)。 这 种 轴 对 称 Stackel 


势 的 第 三 运动 积分 ( 见 文献 图 的 中 表 1 和 文献 回 中 的 式 (16)) 为 : 
(+7)RA LAHE) 1 z ,5 ae | r? 2, Tare 
LN (y - o)(4 — v) 2y-aV up (so ): cas 
XE, (r,z,0) HIER, zo = Ry — o 为 z 轴 上 的 焦距 。 
恒星 盘 分 布 函数 假 定 为 一 个 三 运动 积分 模型 ; 
(8X) E — Ea B 
f Ls) 2 ep (a) x ea 


R3(E, 7 2 d a 
vr ER) G3 
c Bienaymé 等 人 也 是 用 X? 方法 拟 合 参数 。 由 式 (OI KE (从 分 布 函数 得 到 的 速 
o 度 的 若干 阶 平均 )， 在 观测 数据 上 ， 通 过 拟 合 从 三 个 金属 样本 中 的 观测 得 到 矩 ， 即 数 密度 
v(z) 和 垂 向 速度 弥散 02. (z) 来 确定 垂 向 力 的 参数 (只 考虑 了 F, AET RAWX z AN 
的 关联 )。 即 通过 使 用 x? 这 种 积分 数据 处 理 方法 使 建 模 的 量 和 观测 的 量 之 间 的 差异 最 小 化 : 


X! — X2 + x2 t q(V-(9.5 kpc) — V.(7.5kpe))? . (114) 


Bienaymé 4& A.P" 得 到 的 太阳 附近 盘 组 分 的 面 密度 Fuise(Ro) = (44.4 + 4.1) Mo : pc-2， 
太阳 附近 的 暗物质 密度 ppm(z = 0) = (0.0143 + 0.0011) Mo - pc-3，Oort 极限 下 重子 
密度 和 总 密度 分 别 为 pparyons(z = 0) = (0.077 + 0.007) Mo - pc? 和 pparyons(z = 0) = 
(0.091 + 0.005 6) Mo - pc-3。 
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JE Yb, Bienaymé 等 人 ”将 一 般 轴 对 称 引力 势 近似 并 拼凑 为 Stickel 势 形式 ， 即 令 
V(R, z) = V (A, v) 为 在 位 置 (Ri, 21) = (33,14) 处 局 部 近似 的 真实 引力 势 ， 其 Stëkel 势 为 


D(A, v) ( 详 见 其 附录 A)。 并 以 和 前 述 同样 的 方法 基于 Stackel 势 假定 了 I 的 表达 式 : 
1-1 1, 
— V(R,z) E 9": : (115) 


KH, W(R,z)=- Ivan, z) - V(VÀ, 0) ata, 


Bienaymé"" 还 基于 Stäckel 势 近 似 给 出 了 一 般 三 轴 引 力 势 下 的 估算 第 二 、 第 三 运动 积分 
D, I 的 一 种 近似 的 解析 表达 式 ( 见 本 文 区 5 了 小节 )。 然 后 讨论 了 基于 此 的 分 布 函数 模型 : 


f=exp{—n [e+ (bh en) } (116) 


并 用 Jeans 方程 进行 了 检验 。 
322 ”一 个 更 多 参数 的 分 布 函 数 模型 

下 面 介绍 一 种 参数 更 多 、 复 杂 的 分 布 函数 模型 ”™“。Famaey 56 AP" 构建 基于 3 个 运动 
积分 的 分 布 函数 F (E, L,Is)， 从 而 实际 星系 的 恒星 盘 通 过 这 种 形式 的 分 布 函数 (或 者 若干 
此 种 分 布 函数 的 线性 组 合 ) 来 正确 地 拟 合 出 。 他 们 的 方法 是 ，(1) 对 于 DS 的 计算 ， 简 单 套用 
Stickel HAY I. 表达 式 ; (2) 对 于 分 布 函数 的 形式 ， 在 前 人 的 一 个 基于 2 个 运动 积分 的 分 布 
函数 "的 基础 上 ， 引 入 更 多 自由 参数 作 拓展 。 


分 布 函 数 的 具体 形式 为 : 
u E — S,,(Lz) " l 5 
FUE Lu 43 TU) (zz = E) (pt+qE+rL? + sI)? , (117) 
BH, fU.) = PLS. G) Gs s) 909, E, Ley 15 28 3 个 运动 积分 ， 
pq r, s 为 参数 。 根据 这 种 分 布 函数 ， 在 (E, L2) FE, AR z = zo AE, TPR 
频率 ， 则 得 到 边 BAR E >P- Fa (L, =v- 6a): 可 以 用 E= S,,(L,) 表示 包 络 线 。 


通过 计算 分 布 函数 下 若干 阶 速 度 的 平均 ， 即 可 得 到 该 分 布 函数 (运动 积分 空间 ) AFB: 
sn n (A, V) =f f fr (E, Lz, i3) yvy v, dvidv,dv, , (118) 


从 而 0 阶 矩 为 质量 密度 如 p(A, v) = Hooo(A,z)， 二 阶 矩 为 速度 弥散 如 po2 (A, v) = poo, 3(À, v)e 
文献 EA 的 第 4 章 给 出 了 一 种 方法 ， 通 过 由 积分 区 域 边界 的 线性 组 合 表 示 积 分 区 域 中 的 点 ， 
来 计算 其 各 阶 矩 中 的 积分 ， 本 文 附录 B 简介 了 其 基本 思想 。 有 了 和 矩 的 计算 ， 文 献 EB 的 第 5 
章 介 绍 并 分 析 了 分 布 函数 中 各 种 参数 所 体现 的 物理 特征 。 

文献 EQ) 最 后 介绍 了 对 实际 的 恒星 盘 建 模 的 方法 : 通过 将 分 布 函数 组 分 作为 基础 函数 
线性 组 合 。 这 需要 在 配置 空间 中 引入 一 个 网 格 (wi;, zi)， 并 最 小 化 方差 (如 对 于 分 量 Fa 代表 
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0 阶 矩 ， 即 密度 ul) 9(@i, z) 时 ): 


2 
<= >; (>， CALS. Qo (i; 21)) — pal@z, zi) ; (119) 
i A 
其 中 ， A= (a1, Q2, B, Ò, 7T, 70, 8, a) 为 分 布 函数 的 参数 集 ， CA 为 系数 。 拟 合 程序 中 不 断 加 
入 满足 要 求 特征 的 新 组 分 ， 直 至 x? 的 最 小 值 不 显著 变化 。 


4 讨论 和 总 结 


前 面 内 容 聚 焦 于 准 稳 态 时 星系 的 分 布 函数 特点 及 建 模 。 但 是 ， 与 常见 的 稳 态 、 准 稳 态 的 
热力 学 系统 相 比 ， 自 引力 系统 有 些 特点 本 质 不 同 ， 这 些 特点 也 许 会 限制 前 述 方法 的 适用 范 
围 。 另 一 方面 ， 当 考虑 处 于 非 稳 态 (远离 平衡 时 ) 的 星系 结构 和 演化 过 程 ， 或 者 考虑 少 体 引 
JRA (粒子 数 较 少 的 引力 系统 ) 中 的 具体 恒星 的 运动 ， 上 述 基 于 作用 量 的 建 模 方法 是 否 
效 ， 或 者 说 如 何 使 得 上 述 方法 仍 能 应 用 于 某 些 具体 场合 ， 这 是 非常 值得 探索 的 课题 。 

在 本 章 里 ， 我 们 将 首先 在 动力 系统 理论 的 视角 下 (dynamical-systems theory)， 研 究 自 
引力 系统 的 动力 学 ( 见 4.1 节 ); 然后 介绍 和 评述 近年 来 研究 者 把 “作用 量 - 角 变 量 ”方法 应 
到 完全 不 可 积 的 星系 现象 和 星系 结构 的 进展 (共振 和 星系 棒 ， 见 4.2 节 )， 应 用 于 非 平衡 的 
星系 过 程 一 一 潮汐 星 流 (tidal streams) 的 研究 ， 包 括 星 流 的 形成 过 程 ， 在 作用 量 空 间 、 角 变 
量 空间 对 星 流 的 刻画 ， 以 及 把 星 流 用 来 拟 合 主星 系 的 引力 势 分 布 ( 见 4.3 节 ); 最 后 是 全 文 的 
总 结 (JL 4.4 节 )。 

4.1 ”动力 系统 理论 视角 下 的 研究 

从 统计 物理 角度 ， 引 力 N 体系 统 是 一 种 长 程 相互 作用 系统 。 与 一 般 的 近 程 作用 系统 
不 同 ， 这 类 系统 中 任何 粒子 之 间 的 作用 几乎 都 不 能 忽略 (粒子 间 的 作用 势 x r^. Ha 
小 于 位 形 空间 维 数 d;， 假 设 粒 子 在 全 空间 均匀 分 布 ， 那 么 总 作用 势 V x roye 将 是 发 散 
的 )， 并 且 与 势能 有 关 的 宏观 量 不 再 具有 统计 物理 所 经 常 要 求 的 可 加 性 。 另 一 方面 ， 自 引力 
系统 由 于 可 形成 核 - 替 结 构 使 其 人 可 能 不 具有 最 大 值 ， 因 此 并 没有 一 个 真正 的 平衡 态 ;， 一 般 
来 说 ， 相 比 于 近 程 作用 系统 ， 自 引力 系统 的 弛 豫 过 程 可 以 更 快 ， 例 如 剧烈 弛 鸳 机 制 (violent 
relaxation)”™”。 这 些 特点 导致 自 引力 系统 与 统计 物理 中 常见 的 近 程 作用 系统 显著 不 同 ， 往 
往 不 能 用 平衡 态 统计 方法 来 很 好 地 研究 。 

从 动力 学 角度 ， 作 用 量 - 角 变 量 坐标 将 n 自由 度 可 积 系统 的 2n. 维 相 空间 变 成 n 维 环 面 
T"， 如 果 存 在 一 组 非 零 整数 ki, ko, kae, kn WAPA N 满足 y kiQ = 0 时， 则 存 
在 共振 。 共 振 关 系 可 有 不 同 多 种 。 庞 加 莱 (Poincaré) 截面 法 就 是 在 相 空 间 ( 即 坐 标 维 十 速度 
维 ) 中 取 一 个 低 于 相 空间 维度 的 “ 面 >， 记 录 系 统 的 轨道 穿 过 这 个 面 的 点 ， 并 研究 这 些 落 点 
规律 ， 该 方法 是 降 维 研究 运动 的 常见 且 直 观 的 方法 。 以 3 自由 度 可 积 系统 为 例 (如 在 Stáckel 
势 形 式 的 星系 中 的 任意 一 个 恒星 轨道 )， 该 系统 在 相 空间 中 呈现 为 一 个 三 维 环 体 (3-torus); 
若 系统 的 其 中 两 个 频率 比 3/022 为 有 理 数 ， 则 在 其 相应 的 二 维 环 面 的 运动 是 周期 的 ， 其 在 
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环 面 上 的 运动 轨迹 为 周期 轨 线 ， 在 庞 加 莱 截 面 上 则 是 一 些 离散 的 不 动 点 ; 若 运动 是 准 周期 
的 ， 则 是 布 满 环 面 的 准 周期 轨 线 ， 落 点 在 截面 上 将 形成 一 条 闭合 曲线 。 

当 可 积 系统 受 不 可 积 微 扰 且 成 为 近 可 积 系统 ( 近 可 积 系 统 是 大 部 分 的 相 空 间 仍 被 规则 
轨道 占据 的 受 扰 可 积 系 统 )，KAM 定理 给 出 了 近 可 积 系 统 中 仍然 稳定 存在 局 部 的 不 变 环 面 
的 条 件 。KAM 定理 是 20 世纪 中 时 Kolmogorov, Arnold 和 Moser 提出 的 重要 定理 ， 后 续 的 
工作 放宽 了 定理 的 条 件 。 对 于 截面 上 的 不 动 点 ， 当 施加 扰动 ， 简 单 来 说 ， 经 过 多 次 运动 后 ， 
XX dii. (不 稳定 ) 不 动 点 将 逐渐 偏离 ， 椭 圆 (稳定 ) 不 动 点 则 将 围绕 不 动 点 形成 闭合 环 面 ， 或 在 
满足 KAM 定理 的 一 定 区 域内 ( 即 KAM 环 面 ) 稳定 存在 ， 或 形成 更 高 阶 的 新 的 两 种 不 动 点 。 
如 此 达 代 ， 将 导致 更 复杂 和 不 可 预知 的 演化 ， 从 而 产生 混沌 。 此 外 ， 对 于 n > 2 自由 度 的 近 
可 积 系 统 ， 相 空间 为 2n 维 ， 等 能 面 (或 称 能 量 流 形 ) 为 (2n — 1) 维 ， 而 nn 维 的 KAM 环 面 不 
再 能 将 等 能 面 分 割 成 闭 区 域 的 集合 ， 即 存在 阿诺德 扩散 现象 (Arnold diffusion). Iu vage 
BUS n> 3 自由 度 的 动力 系统 的 行为 更 为 复杂 ， 在 稳定 的 轨道 之 间 存 在 着 高 度 不 稳定 的 轨 


道 ， 这 样 运动 轨迹 就 可 以 绕 过 KAM 环 面 到 其 他 区 域 ， 仍 会 出 现 混沌 ， 人 们 对 于 阿诺德 扩散 
ps 的 理解 ， 目 前 仍然 非 党 有限。 大 臻 上， 我 们 可 以 说 ， 有 理 环 面 的 破坏 源 于 运动 模式 的 非 线性 
> 共振 ， 非 线性 共振 进一步 导致 混沌 运动 中 。 现 代 的 “动力 系统 ”学 科 本 身 也 在 拓扑 和 分 析 
+ 上 有 了 更 严密 的 定义 和 更 多 的 理论 ， 我 们 期 待 着 这 些 现代 理论 能 在 自 引力 系统 的 动力 学 研 
co 究 中 大 放 异 彩 。 

= 4.2 ”应 用 于 “不 可 积 ”的 星系 现象 


绝 大 多 数 自 引力 系统 (如 星系 ) 都 不 是 可 积 系 统 ， 即 使 数值 计算 ， 由 于 混沌 的 存在 ， 也 
几乎 不 可 能 精确 计算 某 个 具体 天 体 (如 某 个 恒星 ) 实际 的 长 期 运动 。 尽 管 如 此 ， 对 于 处 于 准 
稳 态 的 恒星 集体 ( 即 星系 ) 而 言 ， 前 述 章节 所 介绍 的 Binney 等 人 的 工作 表明 : 恒星 集体 的 统 
计 性 质 ， 例 如 各 种 分 布 ， 是 可 以 利用 “作用 量 - 角 变 量 ” 方 法 来 可 靠 地 获得 的 。 我 们 进一步 
感 兴趣 的 是 ,“ 作 用 量 - 角 变 量 ” 方 法 是 否 可 以 应 用 于 共振 、 混 沌 等 所 谓 “ 不 可 积 的 ”的 恒星 
运动 ， 可 以 应 用 于 哪些 具体 场合 ? 最 近 Binney ™ ”研究 了 把 环 面 映射 方法 应 用 到 棒 旋 星系 
的 共振 轨道 的 情况 (以 银河 系 观 测 数据 为 实例 )。 他 在 文献 ra] 中 ， 分 别 研究 了 在 银河 系 棒 
的 这 种 扰动 强度 下 ,“ 作 用 量 - 角 变 量 ” 工 具 对 于 外 林 德 布 拉 德 共振 、 共 转 共 振 、 内 林 德 布 拉 
德 共振 的 适用 性 。 他 发 现 ， 对 于 银河 系 被 外 林 德 布 拉 德 共振 和 共 转 共振 所 俘获 的 恒星 轨道 ， 
环 面 映 射 法 能 很 好 地 拟 合 出 ;但 是 更 靠近 棒 内 部 的 内 林 德 布 拉 德 共振 俘获 的 轨道 ， 环 面 映射 
法 就 显得 不 够 好 了 。Binney” 更 进一步 把 上 述 方法 应 用 到 太阳 邻 域 恒星 的 速度 空间 数据 中 ， 
发 现在 银河 系 的 棒 模 型 下 ， 环 面 映 射 法 所 拟 合 的 共 转 共振 情形 与 观测 数据 一 致 ， 而 外 林 德 布 
拉 德 共振 情形 不 被 数据 所 支持 ， 也 就 是 说 ， 基 于 环 面 映射 法 和 微 扰 理论 的 计算 ， 对 于 棒 旋 星 
系 中 的 共振 俘获 现象 ， 可 以 有 很 好 的 实用 性 了 。 

4.3 ”应 用 于 远离 平衡 的 星系 过 程 : 潮汐 星 流 

当 一 个 星团 或 矮星 系 ( 称 为 前 身 天 体 ) 靠近 一 个 大 星系 ( 称 为 主星 系 ) 时 ， 前 身 天 体 各 个 
不 同位 置 所 受 的 引力 不 同 ， 从 而 发 生 形状 和 速度 的 变化 ， 当 所 处 的 引力 场 梯度 足够 大 时 ， 前 
身 天 体 的 恒星 们 会 被 主星 系 巨 大 潮汐 力 逐 步 拉 扯 ， 或 多 或 少 的 恒星 脱离 前 身 天 体 ， 常 形成 长 
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条 状 或 称 纤维 状 的 结构 ， 此 即 潮汐 星 流 ( 详 见 文 献 3 可 )。 

研究 者 发 现 ， 对 于 星 流 的 这 种 形成 过 程 ， 以 及 对 于 星 流 特性 的 刻画 ， 最 方便 的 方式 是 基 
于 “作用 量 - 角 变量 ”体系 ””。 当 前身 矮星 系 或 星团 的 成 员 恒星 没有 被 潮汐 剥离 ， 随 着 前 
身 天 体 在 主星 系 引 力 势 中 运动 时 ， 在 主星 系 参 考 系 中 ， 这 些 恒 星 的 作用 量 是 不 守恒 的 。 但 
且 某 恒星 被 潮汐 力 和 剥离 前 身 天 体 之 后 ， 就 可 以 忽略 前 身 天 体 以 及 星 流 中 的 其 他 恒星 对 它 的 
引力 作用 ， 认 为 只 是 在 主星 系 的 引力 势 下 运动 。 那 么 ， 自 剥离 时 刻 起 ， 该 恒星 的 作用 量 就 
直 保持 剥离 前 的 数值 ， 不 再 随时 间 变 化 了 。 下 面 我 们 沿用 Eyre 和 Binney” 的 记号 和 表述 方 
式 ， 来 简介 这 个 方法 。 een (小 星系 或 星团 ) 位 于 固定 背景 势 下 的 规则 轨道 上 ， 其 作 
HEX Jo。 背 景 势 具有 哈密 顿 量 H(J) KF J 是 潮汐 星 流 中 的 菜 恒 星 的 作用 量 ， 一 般 与 
Jo 相差 不 远 。 对 于 潮汐 星 流 中 的 茶 个 恒星 ， 其 作用 量 是 常量 ， 角 变量 随时 间 线 性 增加 。 在 
Jo 附近 ， 对 哈密 顿 量 作 泰勒 展开 : 


— 


c 


1 
H = Ho + R: ÒT + 50J D.6J | (120) 


其 中 6J =J — Jo; Q% mu ee RUNDE. Qu, = SF s,s 万 是 哈密 顿 量 的 海 森 矩 阵 
( 即 五 的 三 阶 微分 )，D 则 位 于 前 身 天 体 轨道 附近 的 某 恒星 的 轨道 频率 为 : 


ij 一 ara" 
QAIJI)=NM+D-6F . (121) 


关于 星 流 的 成 员 恒星 ， 它 们 的 作用 量 和 和 角 变 量 都 是 有 着 某 种 分 布 ， 即 某 种 程度 的 弥散 ， 
AJ 和 Ab。 由 前 述 可 知 ， 成 员 恒星 的 作用 量 弥散 AJ 自从 星 流 形成 后 就 不 变 了 ; 但 角 变 
弥散 不 会 保持 在 初始 的 和 60， 而 是 将 随时 间 线 性 增 大 。 角 变量 弥散 的 关系 式 为 


iN 
E 


AO(t) — tAQ + A0; ~tAQ ， (122) 


$ 
+t 
N 
a 
a 
am 


由 于 潮汐 流 形成 后 ，0o 将 远 小 于 正比 于 时 间 的 那 一 项 。 于 是 有 


AQ(J) 2 D: AJ . (123) 


上 式 表明 海 森 矩阵 D 的 作用 是 : 把 恒星 们 在 作用 量 空间 的 分 布 (参见 文献 Bo 的 图 10) 线性 
映射 到 角 变 量 空间 的 分 布 (参见 文献 Ba] 的 图 13). D 是 一 个 对 称 和 矩阵 ， 因 此 其 特征 取决 于 
其 相互 正 交 的 特征 向 量 6&, (n = 1,2, 3)， 和 实 特 征 值 ,,”。 海 森 矩 阵 决定 了 星 流 的 拉 伸 方 
向 ， 即 星 流 沿 着 最 大 特征 方向 拉 长 


A0 


什么 样 的 引力 势 分 布 形式 (必要 条 件 ) 可 以 形成 星 流 等 技术 细节 ， 这 里 不 作 详细 介绍 ， 
读者 可 以 参考 Sanders 的 博士 论文 (第 5 章 和 第 6 章 )”。 特 别 是 该 论文 第 6 章 详 细 介绍 了 如 
何 应 用 星 流 数 据 来 拟 合 限定 主星 系 的 引力 势 分 布 。 
最 后 ， 我 们 要 指出 的 是 : 当前 身 天 体 为 很 小 的 星系 或 者 星团 时 ， 星 流 的 轨迹 可 以 粗略 地 
(1 阶 近似 ) 看 作 是 前 身 天 体 在 主星 系 中 的 轨道 ， 因 此 ， 过 去 研究 者 常 简 单 假定 星 流 为 一 条 轨 
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道 ， 从 而 拟 合 主星 系 引力 势 ， 此 即 文献 中 所 说 的 轨道 拟 合法 (例如 文献 BE 到 )。 但 是 ， 基 于 上 


面 的 分 析 ， 星 流 的 轨迹 并 不 必然 与 前 身 天 体 的 轨道 或 者 任何 其 他 轨道 重 


A 
&( 
图 13); 因此 轨道 拟 合法 会 导致 较 大 的 主星 系 引力 势 偏 差 ( 详 见 文献 [20] 的 第 6 章 )。 


4.4 


变量 
APER 


A 


结束 语 


本 文 主要 目的 是 介绍 星系 动力 学 建 模 方 面 的 新 进展 ， 即 作用 量 计 算 方 法 ， 以 及 近年 来 基 
于 作用 量 (或 运动 积分 ) 的 分 布 函数 建 模 在 银河 系数 据 上 的 应 用 ， 另 外 还 简介 了 “作用 量 - 角 
”方法 在 不 可 积 的 星系 现象 (共振 和 星系 棒 ) 和 非 平衡 过 程 (例如 潮汐 星 流 的 形成 ) 的 应 
j， 以 期 为 中 国 在 “Gaia 二 LAMOST” 黄 金 时 代 的 银河 系 研究 做 些 贡献 。 目的 之 二 ， 厘 清 
上 述 方 法 及 其 应 用 的 理论 基础 和 脉络 。 在 实际 或 理想 的 星系 中 ， 唯 一 的 一 类 可 积 系统 是 具 


备 Stückel 形式 引力 势 的 星系 ; 相应 地 ， 几 乎 所 有 星系 建 模 工作 的 解析 性 基础 ， 都 是 Stückel 
势 理论 。 因 此 ， 我 们 还 在 本 文 第 1 章 花 不 少 篇 幅 从 NN 个 粒子 的 经 典 力学 (哈密 顿 方程 和 刘 维 
尔 方程 ) 出 发 ， 依 次 引入 平均 场 近似 ， 介 绍 无 碰撞 玻 尔 兹 曼 方 程 ， 然 后 在 经 典 力学 的 正则 变 


换 以 及 相 空 间 的 语 境 下 ， 引 入 运动 积分 、 正 则 环 面 坐标 即 作 用 量 和 角 变量 ， 接 着 就 介绍 唯 


的 严格 可 积 的 (理想 的 ) “星系 ”， 即 介绍 Stackel 势 理 论 的 各 种 解析 式 。 介 绍 这 些 来 龙 去 脉 ， 


期 望 我 国学 者 不 仅 能 应 用 Binney 等 国外 专家 发 明 发 展 的 方法 来 分 析 银 河 系数 据 ， 还 能 够 在 


理论 基础 和 方法 上 同样 有 所 创新 。 
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Extensions of Stackel Gravitational Potential Theory and 


Their Application in Galactic Modeling 


GU Jian-yu!?, DONG Xiao-bo!? 


(1. Yunnan Observatories, Chinese Academy of Sciences, Kunming 650216, China; 2. University of 
Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: Stáckel gravitational potentials are the most general class of separable potentials. 
Galaxies with a Stáckel form belong to complete integrable systems, in which the stellar or- 
bits are all regular and their integrals of motion (I) can be expressed analytically. Integrals 
of motion, particularly actions (J') as a special kind of integrals, can simplify the description 
of stellar motions. An important approach to study the kinematics and dynamics of stars 
in galaxies is by means of action space; e.g., using actions as the arguments of distribution 
functions (DFs), f(Jr, Lz, Jz). One can estimate actions for the stellar motions in general 
galaxies by locally approximating the gravitational potentials with Stackel models. Recent 
years, there have been important progresses in extending the use of Stáckel potential theory 
in general galaxies: in the respect of computing actions (or integrals of motion) for general 
galaxies, researchers have developed several fast algorithms (e.g., Stáckel Fudge and Stáckel 
Fitting), and several convergent high-precision algorithms based on torus mapping (e.g., 
generating function from orbit integration method and iterative torus construction method), 
as well as formalisms of analytically expressing “approximate integrals of motion" of gener- 
al galaxies (e.g., the quasi-/3); in the respect of constructing DFs for galaxies and galactic 
components, with the actions estimated by the above methods or directly using integrals for- 
mulae of Stáckel potentials (particularly the [3 expression), some DF models are proposed in 
the form of f(Jr, Lz, J2) or f(E, Lz, I3). These theoretical and methodological developments 
make possible DF modelling for galaxies. This article reviews the above developments, and 
introduces several recent representative studies that use action- or integral-based DF models 


to fit the MW data, as well as the recent developments in extending the use of the (action, 
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angle) methodology to non-integrable phenomena (such as galactic resonances and bars) and 
out-of-equilibrium processes (such as the formation of tidal streams), with a hope of being 


helpful for Chinese astronomical community to welcome this “Gaia + LAMOST" era. 


Key words: Stáckel gravitational potentials; action-angle variables estimation; Milky Way 


distribution function modeling 


附录 AO 计算 作用 量 、 角 频率 和 角 变 量 的 步骤 


| 


这 里 介绍 Sanders 和 Binney 的 程序 中 数值 计算 作用 量 、 角 频率 和 角 变 量 的 步骤 (文献 
四 中 第 8 节 及 附录 ; 文献 [22] 的 附录 A)。 程 序 中 的 这 些 中 间 量 的 计算 以 及 轨道 分 类 ， 来 自 


Sanders 和 Binney 的 TACT 程序 ” 的 /tact/aa/inc/stackel_aa.h 中 的 函数 定义 。 
(1) 作用 量 J 
作用 量 的 计算 公式 为 
1 OW, 1 2 T4 
dec = rs dr -— = f pdr = a [p-|dr y TAA, LL, V, (125) 


KP, r, 下 是 友 =0 时 的 根 。 注 意 对 于 有 些 环 轨 道 ，r 的 积分 上 下 限 可 能 是 —o, —8,—y 
具体 可 见 文献 加 引 的 表 1， 其 详细 讨论 可 见 文献 四 中 的 轨道 分 类 内 容 。 
(2) AME 2 
由 哈密 顿 方程 (2; = 6 —0H/8J.) 以 及 微分 关系 ， 可 以 得 到 ”: 


OF _,_ ,0 OS, OJ, 
aE 1 A 3E £2 JE (2, AE 
OE OJ oJ ðJ 
—_ = = Q H Q, v 
ôl "ewm OI * Ol 8L (126) 
OE OJ oJ ðJ 
= = (2 Q H Q v 
Ola 0 入 Ola H Ofs v als 
将 这 个 关于 角 频 率 的 线性 方程 反 解 可 得 角 频 率 ， 如 
TO d) 8(J, Jus Jv) 
= Fath, i) P= XELOLY 127 
* TOES) ODE LIES ， (127) 
可 用 式 (E23) 与 动量 和 运动 积分 的 关系 计算 
OJ, | PUR. I=E, h, 五. (128) 


OI 2n OI 
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HP, H Stackel 势 中 动量 和 运动 积分 的 关系 可 得 : 


Opr »- 

OE l- = Ap. (T +a)(7 + 8)(7 +7) 

Pr ls zm z Opr le (129) 
0j TOE 

Op, 1 Op. 


注意 上 式 右 端 作为 分 母 的 p;， 在 积分 边界 会 消失 。Sanders 和 Binney T FRAM: 


1 1 
T —fsinÜ-T, T= 5(7- +74), 7 = 5(7- T+) ; (130) 


使 被 积 函数 在 积分 边界 更 平滑 地 变 到 零 。 这 样 ， 每 一 点 积分 变量 的 值 都 得 到 了 ， 再 数值 积分 
得 到 角 频 紊 。 程 序 中 用 的 是 高 斯 - 勒 让 德 正 交 格 式 来 进行 积分 。 
(3) AZE 0 
用 哈密 顿 特征 函数 计算 角 变 量 ， 表 示 为 : 
oW OW OI 
um umo (131) 


I: 


计算 OW/OI 所 用 的 W(r, J) A: 
w= SW, 2» op Fo) f dr7’|p’| , (132) 
其 中 ， 被 积 函 数 的 处 理 类 似 于 前 述 所 示 ; 函数 F, (p. c) 是 考虑 轨道 类 型 的 特点 ， 为 消除 椭 


球 坐标 中 的 简 并 性 以 使 角 变 量 履 盖 笛 卡尔 坐标 中 一 个 振动 的 全 部 范围 而 引入 的 (文献 [22] 附 
录 中 的 式 (A6))， 即 计算 中 还 要 进行 轨道 方面 的 判断 和 修正 。 


附录 B  Famaey 等 人 ”的 运动 积分 分 布 函 数 模型 的 矩 的 计算 处 理 方 
法 简介 


对 于 式 (I3)， 从 被 积 的 速度 空间 转换 到 3 个 运动 积分 所 构成 的 空间 


gti f(L E S 
= ) mar? EdI » 
Hl,m,n )* sen L2 Ug d / fa d 3 TJ laze = 3] (133) 


VN RE. (Ix — hy "Ue ; 


# HGB) = A+A) - E] - 55 3. (EL) = (e 9166) - E] - 


N 
4 


202306.00347v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


186 天 文学 进展 39 3$ 


YEY 72, dE (E, I3) 面 的 积分 区 域 为 : 


I; (BE,L:) < Is < I3 (E, L2) 


S,(L,)&EX (134) 


p+qE+rL?+sIz 20 


通过 对 积分 边界 线 (上 式 取 等 号 时 ) 的 线性 组 合 ， 来 表示 式 (OR) 中 的 积分 区 域 中 的 点 ， 
E -— Sa 
So m Szo 
之 后 ， 将 (EB, 13) 的 积分 区 间 转 换 到 新 的 积分 变量 空间 (固定 工 :)， 然 后 对 比 变换 前 后 的 每 一 
项 ， 解 出 上 式 中 的 组 合 系数 v, u, 二 的 关系 ， 从 而 用 一 维 数值 积分 计算 其 矩 (更 多 细节 见 文献 
[22] 中 的 4.2 和 4.3 节 )。 


= t(j — I) + ull — I5) - v(p--qE  rL2 + sI). (135) 


